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Introduction 


The NASA Langley Research Center 0.3-m Transonic Cryogenic Tunnel (TCT) has 
been operational since 1973. It is a closed circuit tan driven pressure tunnel in which the heat 
generated by the fan operation is cancelled by injection of liquid nitrogen (LN 2 ) into the 
tunnel. The injected LN 2 evaporates into the tunnel and cools the tunnel resident nitrogen gas 
(GN,) mass. The gaseous mass buildup and pressure variation caused by LN Z injection is 
corrected by controlled discharge of warmer GN 2 tunnel gases. The gas temperature, pressure, 
and flow Mach number of the tunnel can be controlled by adjusting the LN 2 injection and the 
GN Z discharge for a given fan speed. 

In 1988 a microcomputer was installed at the 0.3-m TCT to control the temperature, 
pressure, LN, back pressure, fan speed, Mach number and Reynolds number. This system 
replaced the existing microprocessor control system. This new controller system consists of 
a PC AT clone with a 12 MHz CPU, an EGA video, a hard disk, a floppy disk, and an 8 
channel digital to analog converter. This Controller was designed specifically for ease in 
modification and operational use. This system has now been operating tor several thousand 
hours providing a sate anil ellicient control ot the tunnel. (Ret. 1 and 2) 

Since the Controller has been in use, the staff has gained an extensive knowledge of 
the Controller and its features. Through the operational use of the Controller several 
recommendations for improvements have been suggested. Many of these suggestions have 
been incorporated into the Controller and are described later in this document. The control 
laws and operations of the Controller have not changed from the original Controller. 

The purpose of this document is to provide the source code listing of the 0.3-m TCT 
Controller which is in Appendix A. Also included is the source code listing for the tunnel 
simulator (Appendix B) and a complete listing of the variables used in each source code 

(Appendix C). 
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Controller Changes 


Described below are the changes made to the original Controller. A reason and a 
description of the change are given for a better understanding of the improvements made to 
the Controller. Changes to the Controller software should be performed by a person who fully 
understands the BASIC language. This should be performed only with the knowledge and 
approval of the Facility Safety Head of the ().3-m TCT. 

Set Tunnel Operation Limits 

A requested change was to have the ability for the user to set the tunnel operating 
limits. An example of this is to have the user set a maximum pressure limit of 50 psia. This 
limit is below the tunnel maximum operating pressure of 88 psia and is therefore an acceptable 
limit. With the user set 50 psia limit, the user would not be allowed to input pressures above 
50 psia while operating the Controller. This ability to have the user set the tunnel operating 
limits required an addition to the beginning of the original Controller. 

The Controller is now designed so that when the Controller is first executed a 
Maximum - Minimum Limits program is run which allows the user to set the operating limits 
of the tunnel. The user set limits must be within the operating envelope of the tunnel. This 
Maximum - Minimum Limits monitor display is shown in figure 1. This display is self- 
explanatory and operates using similar keystrokes as the Controller. (Ref. 2) The Maximum - 
Minimum Limits program allows the user to set a minimum temperature, maximum pressure, 
maximum fan speed, maximum Mach number, and maximum Reynolds number. The user can 
either start the Controller with the tunnel limits or the user set limits. 

To set the pressure limit the user would strike the "P" key on the keyboard. This 
would cause a ",psia" to appear on the screen in the Set Maximum column. The user keys 
in the desired pressure limit and strikes the "Enter" key. If the desired limit is unusable, the 
program will warn the user of the problem and will not allow the Controller to be started. 
The other limits may be set by striking the appropriate T, N, M, or R keys on the keyboard. 
The Maximum - Minimum Limits program has several safety features which prevent the user 
from setting unacceptable or unsafe limits. 

These user set limits are only enforced when the tunnel is in the automatic control 
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mode of the variable. For example, if the user set limit for pressure is 50 psia and is 
operating in automatic pressure control, the tunnel Controller will not allow pressure inputs 
above 50 psia. If the user is operating in automatic Reynolds number control, the user again 
will not be able to exceed the set Reynolds number limit. If the user is operating in manual 
GN Z exhaust valve control, the pressure can exceed the 50 psia limit but not the tunnel 
pressure limit of 88 psia. 

Because the Reynolds number control is obtained through pressure control, the Reynolds 
number will be restricted by the pressure limit of 50 psia and, therefore limits the Reynolds 
number even though a user set limit for Reynolds number is not set. This Reynolds number 
limitation can be easily overcome by lowering the tunnel gas temperature. 

Display of Input Limits 

With the ability of the user to set the tunnel operating limits the user needs to know 
these limits. These limits are now displayed in the bottom two rows of the Controller display 
screen in the section titled LIMITS as shown in figure 2. The limit ranges for each variable 
are displayed within their corresponding control loop columns. 

Improve Screen Layout 

The Controller display layout shown in figure 2 is nearly identical to the original 
Controller except for LIMITS displayed at the bottom two rows of the screen. Improvements 
have been made in color and spacing for the visual warnings for an "Emergency Stop". Also, 
displays are now centered in their perspective control loop columns. 

The Reynolds number display has been changed to allow Reynolds numbers of over 
99 million. Before, when the Reynolds number reached 100 million, the display showed a % 
symbol indicating a display overflow. 

The RPM/MACII LOOP now clears the old set point when changing from manual to 
automatic control. The Pt/Re LOOP also clears the old set point when changing from 
automatic pressure control to automatic Reynolds number control. 
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Aerodynamic Chord Inpul 

In the original Controller chord length input was in meters. This has been changed so 
that the chord input is now in inches. This change has been made to improve ease of use for 
the Controller. 

Tunnel Operating Envelope Safety Catches 

Included in the Controller are safety catches which will safely stop the tunnel by an 
"Emergency Stop" if the tunnel limits for pressure or normal fan speed are exceeded. In the 
original Controller it was possible to over pressurize the tunnel when operating with manual 
GN Z exhaust valve control. These safety catches are additional safety features which are 
duplicated in several other tunnel safety devices separate from the Controller. 

Identify Causes of an Emergency Stop 

The Controller "Emergency Stop" procedure brings the tunnel to a safe condition. The 
original Controller displayed only "Emergency Stop" and if the problem was a sensor failure, 
"Sensor Failure" was also displayed. This was frustrating for the user not knowing what had 
caused the "Emergency Stop". This "Emergency Stop" condition now displays the cause of 
the "Emergency Stop" and identifies the failed sensor during a "Sensor Failure". Causes of 
an "Emergency Stop" are displayed as "Pressure Limit", "Fan RPM Limit", or Sensor 
Failure". During a "Sensor Failure" the failed sensors are identified as "Gas Temperature", 
"Wall Temperature", "Total Pressure", "Static Pressure", and "Screen Pressure". These displays 
will allow the ().3-m TCT staff to correct tunnel problems quickly. 

Improved Organization of Source Code 

The source code has now been divided into logical sections which pertain to certain 
functions of the Controller. These divisions are designed to allow for an easier understanding 
of the source code. Some sections have been rearranged to reduce the number of lines and 
improve the Controller execution times. 
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Improved Cooldown Efficiency 

The automatic cooldown of the Controller is one of its best features. It allows for a 
safe, automatic, and efficient cooldown of the tunnel to any desired temperature set point. 
The temperature control loop achieves this cooldown by using a "use 1 temperature set point 
which is based on an acceptable temperature difference between the tunnel wall and gas. 

During tunnel operations the temperature set point is changed frequently by the user 
and the Controller always uses the previous cycle "use" temperature set point as a starting 
point to reach the desired temperature set point. This works well until the tunnel is already 
*Cooi and a' 'user chooses a warmer' "temperature set point and then changes back to a cold 
temperature set point. Because the Controller uses the previous cycle use temperature set 
point to start cooling, it can lake several cycles before the use tempeialure set point reaches 
the actual tunnel gas temperature. This provides a period during a cooldown in which the 
tunnel is actually warming up waiting for the "use" temperature set point to reach the actual 
tunnel gas temperature. 

This problem has been corrected by setting the "use" temperature set point previous 
cycle value, STP, to the tunnel gas temperature, IT. This occurs only during a new 
temperature input to cooldown the tunnel. 

Improved "On" Set Point Flags 

The "On" set point flags are black bands which are shown when the tunnel is within 
a given tolerance of a desired set point. These "On" set point flags work well to indicate 
quickly and easily when the tunnel had reached a desired condition. However, there was 
confusion about the temperature and Reynolds number "On" set point Hags. 

In the original Controller the temperature "On” set point flag was displayed when the 
absolute difference between the gas temperature, IT, and the "use" temperature set point, ST, 
was less then or equal to 0.3. This condition for displaying the "On" set point flag allowed 
the flag to be displayed when the tunnel was not near the desired temperature set point, ST1. 
This caused some confusion in determining if the tunnel was actually "On" temperature set 
point. To remove this confusion the "On" set point flag for temperature is now displayed 
when the absolute difference between the tunnel temperature, TT, and the desired temperature 
set point, ST1 is less than or equal to 0.3. 
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When operating in the Reynolds number control loop the user chooses a desired 
Reynolds number and a pressure is calculated to give this desired Reynolds number. The 
confusion occurs because Reynolds number control is still through the pressure control loop. 
The "On" set point flag in the original Controller is for pressure only. This can cause the 
"On" set point flag to be "On" when in the Reynolds number control loop while there is still 
a large difference between the desired Reynolds number set point, SRE, and the tunnel 
Reynolds number, RE. The solution to eliminate this confusion was to make a separate 
condition for the "On" set point flag for use when in the Reynolds number control loop. This 
has been accomplished by the addition of the logic flag FL7. The Reynolds number "On" set 
point flag is now displayed when the difference between the desired Reynolds number set 
point, SRE, and the tunnel Reynolds number, RE, is less than or equal to 0.05. 

The controller now has "On" set point flags which display when the conditions of the 
tunnel are those of the desired user set points. 

Rewrite Source Code 

The original source code was written in the BASIC language which requires a line 
number for every line of code used in the program. This line numbering causes some 
difficulty when trying to add to the program. The source code has been rewritten and 
compiled using MicroSoft QuickBASIC (Ref. 3). The advantages of using QuickBASIC are 
that it does not require line numbers for each line of code and all statements are identical to 
the original source code. Line numbers in QuickBASIC are only required on lines which are 
called by other lines. This lack of line numbers makes additions or changes much simpler. 
QuickBASIC also has a BASIC compiler that provides a faster compiled Controller program. 
The source code listings in Appendix A and B are in QuickBASIC. 
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Figure 2. Typical Controller Display. 
















Maximum - Minimum Limits for 0.3-m TCT Controller 
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Figure 1. Maximum - Minimum Limits Display. 
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Appendix A 

Controller Source Code Listing 


’0.3-m Transonic Cryogenic T unnel Controller. 
CLS 
CLEAR 


’Set up of the Digital to Analog Converter. 

DEFINT'I, Z 
DIM E(8), DAC(8) 

ZASE.ADDRESS = &1I2EC 

ZOMMAND.REGISTER = ZASE.ADDRESS + 1 

ZTATUS. REG I STER = ZASE.ADDRESS + 1 

ZATA.REGISTER = ZASE.ADDRESS 

ZOMMAND.WAIT = &1I4 

WRITE. WAIT = &II2 

READ. WAIT = &II5 

ZERROR = &II2 

ZCLEAR = &H1 

ZADIN = &HC 

ZSTOP = &HF 


’Digital to Analog Conversion Board Test. 

OUT &II225, 0 
EOR JJ = 1 TO 300 
BSD = SQR(7) 

NEXT JJ 

AEMP = INP(&II225) 

IE AEMP = 4 THEN 2 

PRINT "Digital to Analog Conversion Problem!": PRINT 
GOTO 10000 

2 OUT &I1224, &II6 

’Program Constants. 

DEL = .1 
CGV = 8 
CII = .18 

CHIN = .18 * .0245 
CLQV = 4 
IW = 39 
K = 1 
KDL = 0! 

KDP = 0! 

KDT = 0! 

KIL = .02 
KIM = 0! 

KIN = A 
KIP = .05 
Krr = .1 
KMM = .3 
KP = 3 
KPL = .2 
KPM = 4.5 
KPN = .6 
KPP = 1 ! 

KPT = 1 ! 
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KT = .04 
MXT = 40 
SLQSC = 17! 

STP = 300! 

XD1.P = 1! 
XIIU'M = 1280 
XPLQ = 5.103429 
XPP = 1.366 
XPS = 1.366 


’Maximum - Minimum Limits Screen I-ayoul. 

COLOR 14, 9 
CLS 

IOC Alii 3, 16: PRINT " Maximum - Minimum Limits lor 0.3-m 1(71’ Controller 
COLOR 12, 9 

LOCATL 5, 25: PRINT 'Tunnel" 

LOCATL 6, 24: PRINT "Maximun" 

LOCATL 5, 38: PRINT “Set" 

LOCATL 6, 36: PRINT "Maximun" 

LOCAITi 5, 60: PRIN T "Set" 

LOCAT L 6, 58: PRIN T "Minimum" 

LOCATL 5, 72: PRINT "Tunnel" 

LOCATE 6, 71: PRIN T "Minimum" 


23, 


COLOR 7, 9 

LOCATE 7, 3: PRINT Temperature" 

LOCATE 8, 4: PRIN T "(kelvin)" 

LOCATE 10, 3: PRIN T "Pressure” 

LOCATE 11, 4: PRIN T "(psia)" 

LOCATE 13, 3: PRIN T "Normal Tan Speed” 

LOCATE 14, 4: PRINT "(rpm)" 

LOCATE 16, 3: PRINT "Mach Number" 

LOCATE 19, 3: PRIN T "Reynolds Number" 

LOCATE 20, 4: PRINT - "(Miln)" ... , .. „ 

LOCATE 22, 3: PRINT M Acccpl runnel Max./Min. Limits and Start 0.3-m ICi Controller. 

3: PRINT "Accept Set Max./Min. Limits and Start 0.3-m TCI Controller." 


COLOR 15, 9 

LOCATE 7, 51: PRINT "T’ 

LOCATE 7, 49: PRIN T C1IR$(24.3) 
LOCATE 7, 53: PRIN T C11R$(243) 
LOCATE 10, 51: PRIN T "P" 

1 OCATE 10, 49: PRIN T (TIR$(243) 
LOCATE 10, 53: PRIN T CNR$(243) 
LOCATE 13, 51: PRIN T "N" 
LOCAITI 13, 49: PRIN T CTIR$(243) 
LOCATL 13, 53: PRIN T C1IR$(243) 
LOCATE 16, 51: PRINT ”M" 
LOCATE 16, 49: PRINT CIIR$(243) 
LOCATE 16, 53: PRINT - CIIR$(243) 
LOCATE 19, 51: PRINT - "R" 
LOCATE 19, 49: PRINT CIIR$(243) 
LOCATE 19, 53: PRINT C11R$(243) 

COLOR 14, 9 
LOCATE 7, 3: PRINT "T” 


IjOCAIE 
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LOCA'Ili 10, 3: PRINT "P" 
LOCA1T 13, 3: PRINT "N" 
LOCATE 16, 3: PRINT "M" 
LOCATE 19, 3: PRINT "R" 
LOCATE 22, 3: PRIN’ T "A" 
LOCATE 23, 10: PRINT "S" 


\innel Maximum - Minimum Limits. 
MAXTl = 340! 

MINT1 = 77! 

MAXP1 = 88! 

MINP1 = 14.696 
MAXN1 = 5600! 

MINN1 = 0! 

MAXM1 = .995 
MINM1 = .15 
MAX RE 1 = 50! 

MINRE1 = 1! 

MAXC11IN = 15.75 
MINCIIIN = .394 
MAXCII = .4 
MINCII = .01 
MAXLQSC = 150! 

MAXSPR = MAX PI / 14.696 
MINSPR = MINP1 / 1 4.696 


MAXT = MAXTl 
MINT = MINI! 

MAXP = MAXIM 
MINP = MINIM 
MAXN = MAXN1 
MINN = MINN1 
MAXM = MAXM1 
MINM = MINM1 
MAX RE = MAXRE1 
MINRE = MINRE1 

COLOR 10, 9 

LOCATE 7, 25: PRINT USING "tf##.#"; 
LOCATE 7, 37: PRINT USING "###.#"; 
LOCATE 10, 25: PRINT USING 
LOCATE 10, 59: PRINT USING 
LOCATE 13, 25: PRINT USING "####." 
LOCATE 13, 60: PRINT USING "####." 
LOCATE 16, 25: PRINT USING "#.###' 
LOCATE 16, 60: PRINT USING "#.###" 
LOCATE 19, 25: PRINT USING "##.##" 
LOCATE 19, 60: PRINT USING 
LOCATE 7, 69: PRINT "Saturation" 
LOCATE 8, 69: PRINT "Temperature" 
LOCATE 10, 71: PRINT "Ambient" 
LOCATE 11, 71: PRINT "Pressure" 
LOCATE 13, 73: PRINT USING "####." 
LOCATE 16, 73: PRINT USING "#.###" 
LOCATE 19, 73: PRINT USING "##.##" 


MAXTl 
MAXTl 
MAXPI 
; MINIM 
; MAXN1 
; MINN! 

; MAXM I 
; MINM1 
; MAXREI 
; MINREl 


; MINN1 
; MINM1 
; MINREl 
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’Maximum - Minimum Limits Inputs. 


A$ = INKEYS 
IF AS 
IF AS 
IF AS 
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IF 

GOTO 5 


AS = 
AS = 
AS = 
AS = 


"T" OR AS = "t" TULN GOSUB 10 
"P" OR AS = "p" THEN GOSUB 100 
"N" OR AS = "n" THEN GOSUB 200 
"M" OR AS = "m” THEN GOSUB 300 
"R” OR AS = ”r" 'I'll LIN GOSUB 4(X) 
"A" OR AS = "a" THEN 500 
"S" OR AS = "s" THEN 600 


’1’cmpcralurc Minimum Limit Subroutine. 

10 COLOR 15, 9 
BEEP 
LR = 7 
LC = 58 

LOCATE LR, LC + 1: PRINT " ,K" 


20 

30 

40 


XX = 0 
B$ = INKEYS 

IF B$ = "D" OR B$ = 
LOCATE LR, LC + 1: PRINT 
RETURN 

"• B$ = "0" 


"d” THEN 30 El ,SE 40 


OR B$ = "1" OR BS = M 2" OR B$ = "3'' OR B$ = "4" OR B$ = "5" OR BS = 


"6" OR BS = "7" OR BS = "8" OR B$ = "9" OR BS = THEN 50 ELSE 20 


50 


60 


70 

80 


XX = XX 
II 


+ 1 
XX 
IF XX 
IF XX 
IF XX 
IF XX 
IF XX 


THEN 111 = VAliBS): LOCATE LR, LC + 1: PRINT USING III 

THEN 112 = VAI/BS : LOCATE LR, LC + 2: PRINT USING H2 

THEN 1 13 = VAI/BS): LOCATE LR, LC + 3: PRINT USING 

THEN LOCATE LR, LC + 4: PRINT CIIR$(46) 

5 THEN 115 = VAI/BS): LOCATE LR, LC + 5: PRINT USING 

5 'HI LIN 60 ELSE. 20 
r 112 * 10 + 113 + 1 15 / 10 


1 13 
115 


AA = III * 100 
MINT = AA ... 

IF MINT >= MAXT1 THEN MINT = MAXT1 
IF MINT <= MINT1 THEN MINT = Ml NTT 
C$ = INKEYS 

IF C$ = CIIR$(13) THEN 80 ELSE 70 

BEEP 

COLOR 14, 9 

LOCATE LR, LC + 1: PRINT USING "###.# MIN'1 
RETURN 


’Pressure Maximum Limit Subroutine. 
KM) COLOR 15, 9 
BEEP 
LC = 36 


110 

120 


LR = 10 

LOCATE LR, LC + 1: PRIN T " ,psia" 
XX = 0 


13$ = INKEYS 

IF BS = "D” OR BS = "d" THEN 120 EINE 130 
LOCATE LR, LC + 1: PRINT " " 

RETURN 


130 

140 


IF BS = "0" OR BS = "1" OR BS = "2" OR B$ = "3" OR BS = "4" OR BS - "5” OR BS 
"6" OR BS = "7" OR BS = "8" OR B$ = "9" OR BS = H . M THEN 140 EI^E 110 

XX = XX + 1 


13 


IF XX = 1 THEN III = VAIXB$): LOCATE LR, LC + 1: PRINT USING III 

IF XX = 2 THEN 112 = VAI<B$): LOCATE LR, LC + 2: PRINT USING 112 

IF XX = 4 THEN 1 13 = VAI.(B$): LOCA TE LR, LC + 4: PRINT USING 1 13 

IF XX = 3 THEN LOCATE LR, LC + 3: PRINT CIIR$(46) 

IF XX = 5 THEN 115 = VAl^BS): LOCATE LR, LC + 5: PRINT USING 115 

II- XX = 5 THEN 150 ELSE 110 
150 AA = 111 * 10 + 112 + 1 13 / 10 + 115 / 1 (X) 

MAXP = AA 

IF MAXP <= MINP1 THEN 160 ECSE 170 
160 MAXP = MINIM 

COLOR 20, 8: LOCATE LR + 1, 36: PRIN T " Warning! Set Max. = Tunnel Min. " 

III) = I 
GOTO 180 

170 COLOR 14, 9: LOCATE LR + 1, 36: PRINT " 

FLO = 0 

180 IF MAXP > MAXIM THEN MAXP = MAXIM 

ion - INKFYIt 

IF C$ = CIIR$(13) THEN 195 ELSE 190 

195 BEEP 

COLOR 14, 9 

LOCATE LR, LC + 1: PRINT USING "##.## MAXP 
RETURN 

’Normal Fan Speed Maximum Limit Subroutine. 

200 COLOR 15, 9 

BEEP 
LC = 36 
LR = 13 

LOCATE LR, LC + 1: PRINT " ,rpm” 

XX = 0 

210 B$ = INKEYS 

IF B$ = "D" OR B$ = "d" THEN 220 ELSE 230 
220 LOCATE LR, LC + 1: PRIN T " ” 

RETURN 

230 IF 13$ = "0" OR 13$ = "1” OR 13$ = "2" OR 13$ = "3" OR B$ = "4" OR 13$ = "5" OR B$ 

’'6" OR 13$ = "7" OR 13$ = "8" OR 13$ = "9” OR 13$ = THEN 240 ELSE 210 
240 XX = XX + 1 

IF XX = 1 THEN III = VAI^BS): LOCATE LR, LC + 1: PRINT USING HI 

IF XX = 2 IlIEN 112 = VAl^BS): LOCATE LR, LC + 2: PRINT USING 112 

IF XX = 3 '11IEN 1 13 = VA1<13$): LOCA TE LR, LC + 3: PRINT USING H3 

IF XX = 5 THEN LOCATE LR, LC + 5: PRINT CIIR$(46) 

IF XX = 4 THEN H5 = VAI^BS): LOCATE LR, LC + 4: PRINT USING H5 

IF XX = 5 IlIEN 250 EI>SE 210 
250 AA = III * 1000 + 112 * 100 + 113 * 10 + 1 15 
MAXN = AA 

IF MAXN <= MINN1 THEN 260 ELSE 270 
260 MAXN = MINN1 

COLOR 20, 8: LOCATE LR + 1, 36: PRINT " Warning! Set Max. = Tunnel Min. " 

FID = 1 
GOTO 280 

270 COLOR 14, 9: LOCA'Ili LR + 1, 36: PRINT " 

FID = 0 

280 IF MAXN > MAXN1 IlIEN MAXN = MAXN1 

290 C$ = INKEY$ 

IF’ C$ = CIIR$(13) IlIEN 295 ELSE 290 

295 BEEP 

COLOR 14, 9 


14 



; MAXN 


LOCATE LR, LC + 1: PRINT USING "####. 
RETURN 


’Mach Number Maximum Limit Subroutine. 

3(K) COLOR 15, 9 
BEEP 
LC = 36 
LR = 16 

LOCATE LR, LC + 1: PRINT " ,Mach” 

XX = 0 

310 B$ = INKEY$ 

IF B$ = "D" OR B$ = "d” THEN 320 FUSE 330 
320 LOCATE LR, LC + 1: PRINT " " 

RETURN 

330 IF B$ = "0" OR B$ = "I" OR B$ = "2" OR B$ = "3" OR B$ = "4" OR BS = "5" OR B$ 

"6" OR B$ = "7" OR B$ = ”8" OR B$ = "9" OR B$ = THEN 340 EKSE 310 
340 XX = XX + 1 

IF XX = 1 THEN III = VAI.(B$): LOCATE LR, LC + 1: PRINT USING III 

IP’ XX = 3 THEN 112 = VAI.(B$): LOCATli LR, LC + 3: PRINT USING 112 

IF XX = 4 THEN 113 = VAI^BS): LOCATE LR, LC + 4: PRINT USING 113 

IF XX = 2 THEN LOCATE LR, LC + 2: PRINT CIIR$(46) 

IF XX = 5 THEN 1 15 = VAGITS): LOCATE LR, LC + 5: PRINT USING H5 

IP' XX = 5 THEN 350 EKSE 310 
350 AA = III + 1 12 / 10 + 113 / 100 + 115 / 1000 
MAXM = AA 

IP' MAXM <= MINM1 THEN 360 ELSE 370 
360 MAXM = MIN Ml 

COLOR 20, 8: IjOCATE LR + I, 36: PRINT " Warning! Set Max. = Tunnel Min. " 

P'IjO = 1 
GOTO 380 

370 COLOR 14, 9: LOCATE LR + 1, 36: PRINT " 

FLO = 0 

380 IF MAXM > MAXM1 IlIEN MAXM = MAXM1 

390 C$ = INKEYS 

IF C$ = CIIR$(13) THEN 395 ELSE 390 

395 BEEP 

COLOR 14, 9 

LOCATE LR, LC + 1: PRINT USING "#.### MAXM 
RETURN 

’Reynolds Number Maximum I.imil Subroutine. 

400 COLOR 15, 9 
BEEP 
LC = 36 
LR = 19 

LOCATE LR, LC + 1: PRINT " ,Miln" 

XX = 0 

410 BS = INKEYS 

IF BS = "D" OR BS = "d" THEN 420 ELSE 430 
420 LOCATE LR, LC + 1: PRINT " " 

RETURN 

430 IP’ BS = "0" OR BS = "I" OR BS = "2" OR BS = "3” OR BS = "4" OR BS = "5" OR BS 

"6" OR BS = "7" OR BS = "8" OR BS = ”9" OR B$ = 'HIEN 440 EKSE 410 
440 XX = XX + 1 

IF XX = 1 THEN III = VAI^BS): LOCATE LR, LC + 1: PRINT USING HI 

IF XX = 2 THEN 112 = VAI^BS): LOCATE LR, LC + 2: PRINT USING H2 
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IF XX = 4 THEN 113 = VAI^HS): LOCATE LR, LC + 4: PRINT USING "#"; 113 
IF XX = 3 THEN LOCATE LR, LC + 3: PRINT CIIR$(46) 

IF XX = 5 THEN 1 15 = VAl^BS): LOCATE LR, LC + 5: PRINT USING H5 
IF XX = 5 THEN 450 EESE 410 
450 AA = III * 10 + 112 + 113 / 10 + 115 / 100 
MAXRE = AA 

IF MAXRE <= MINRE1 THEN 460 EI.SE 470 
460 MAXRE = MINRE1 

COLOR 20, 8: LOCATE LR + 1, 36: PRINT " Warning! Set Max. = Tunnel Min. " 

FIjO = 1 
GOTO 480 

470 COLOR 14, 9: LOCATE LR + 1, 36: PRIN T * 

FLO = 0 

480 IF MAXRE > MAXRE1 THEN MAXRE = MAXRE1 

490 C$ = INKEYS 

IF C$ = CIIR$(13) THEN 495 EliSE 490 

495 BEEP 

COLOR 14, 9 

LOCATE LR, I.C + 1: PRINT USING "MM MAXRE 
RETURN 


’Starts the 0.3-m TCT Controller with Tunnel Maximum - Minimum Limits. 
500 XX = 0 

GO TO 800 


’Starts the 0.3-m TCT (Controller with Set Maximum - Minimum Limits. 
600 IF FIjO - 1 THEN GOTO 5 

XX = 0 

MAXSPR = MAXP / 14.6% 

MINSPR = MINI 5 / 14.6% 

GOTO 800 

’0.3-m 'I'd' T-P/Rc-M Controller Program. 


’Controller Screen I ay out. 

800 COLOR 0, 0 
CI.S 

COLOR 15, 4 

LOCATE 1, 26: PRINT " 0.3-m TUNNEL T-P/R-M CONTROLLER " 


I'OR I = 2 TO 25 
COLOR 15, 1 

LOCATE I, 2: PRINT STR1NG$(10, 0); 

IF' I = 4 11IEN LOCA'IE 4, 2: PRINT STRING$(10, 205); 

IF I = 14 THEN LOCATE 14, 2: PRIN T STRING$(10, 205); 
IF I = 19 THEN LOCATE 19, 2: PRINT STRING$(t0, 205); 
IF I = 23 THEN LOCA'IE 23, 2: PRINT STR1NG$(18, 205); 

COLOR 15, 2 

LOCATE I, 12: PRINT STRING$(12, 0); 

IF I = 4 THEN LOCATE 4, 12: PRINT STRING$(12, 205); 

IF I = 14 THEN LOCATE 14, 12: PRINT STRING$(12, 205); 
IF I = 19 THEN LOCATE 19, 12: PRINT STRING$(12, 205); 
IF I = 23 THEN LOCATE 23, 12: PRINT STRING$(18, 205); 


COLOR 15, 3 

LOCATE I, 24: PRINT STRING$(17, 0); 
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IF I = 4 THEN I DC ATE 4, 24: PRINT STRING$(17, 205); 

IF I = 14 THEN LOCATE 14, 24: PRINT STRING$(17, 205); 

IF I = 19 THEN LOCATE 19, 24: PRINT STRING$(17, 205); 

IF I = 23 THEN LOCATE 23, 24: PRINT STRING$(18, 205); 

COLOR 15, 5 

LOCATE I, 41: PRINT STRING$(18, 0); 

IF 1 = 4 THEN LOCATE 4, 41: PRINT STRING$(18, 205); 

IF I = 14 THEN LOCATE 14, 41: PRINT STRING$(18, 205): 

IF I = 19 THEN LOCATE 19, 41: PRINT SITING $(18, 205); 

IF I = 23 THEN LOCATE 23, 41: PRINT STRING$(18, 205): 


COLOR 15, 6 

LOCATE I, 59: PRINT STR1NG$(18, 0); 

IF I = 4 THEN LOCATE 4, 59: PRINT STRING$(18, 205); 

IF I = 14 THEN LOCATE 14, 59: PRINT STR1NG$(18, 205) 

IF I = 19 THEN LOCATE 19, 59: PRINT STRING$(18, 205) 

IF I = 23 THEN LOCATE 23, 59: PRINT STRING$(18, 205) 

NEXT' I 


COLOR 15, 2 

LOCATE 2, 14: PRINT' " LN PUMP ": PRINT ; 
LOCATE 3, 14: PRINT' " AUTO ": PRINT' ; 

LOCAT E 5, 18: PRINT' ",psi;i": PRINT' ; 

LOCATE 9, 18: PRINT' ",psia": PRINT ; 

LOCATE 12, 18: PRINT' ",%opn": PRINT ; 

LOCATE 15, 12: PRINT' " =": PRINT’ ; 

LOCATE 5, 13: PRINT USING "###.#"; SLQSC 

COLOR 14, 2 

LOCATE 15, 12: PRINT’ "H": PRINT ; 

COLOR 15, 3 

LOCATE 2, 28: PRIN T "TEMP LOOP": PRINT ; 
LOCATE 3, 29: PRIN T " MANUAL ": PRINT' ; 
LOCATE 5, 31: PRINT B ,K(Final)": PRINT ; 

LOCATE 6, 31: PRINT' ",K(Usc) H : PRINT ; 

LOCA TE 9, 31: PRINT ",K-GN2": PRINT ; 

LOCATE 10. 31: PRIN T ".K-WALL": PRINT ; 

LOCATE 12, 31: PRINT ",%opn": PRINT ; 

LOCATE 15, 25: PRIN T " cmp=": PRINT ; 

LOCATE 16, 25: PRINT "A Q%=": PRINT ; 

LOCATE 20, 25: PRIN T "GRAD* ,K/mt": PRINT ; 
LOCATE 21, 25: PRINT' " SAT'= ,K": PRINT' ; 
LOCATE 24, 32: PRIN T "T"; 

LOCATE 24, 30: PRIN T CHR$(243); 

LOCATE 24, 34: PRINT CIIR$(243); 

LOCATE 24, 26: PRINT USING "###"; MINT; 

LOCATE 24, 36: PRIN T USING "###"; MAXI’; 

COLOR 14, 3 

LOCATE 15, 25: PRINT "T": PRINT ; 

LOCATE 16, 26: PRINT "L": PRINT ; 

COLOR 15, 5 

LOCATE 2, 45: PRINT "Pt/Rc LOOP": PRINT ; 
LOCATE 3, 46: PRINT " MANUAL ": PRINT ; 
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LOCATE 5, 48: PRINT ",Psia": PRINT ; 

LOCATE 6, 48: PRINT ”,Miln": PRINT ; 

LOCA'Ili 20, 42: PRINT!' "CIIORD= 7.09,in": PRINT ; 
LOCATE 9, 48: PRINT ”,Psia": PRINT ; 

LOCATE 21, 42: PRINT "P sl= ,psia": PRINT ; 

LOCATE 22, 42: PRINT "Del P= ,psi": PRIN T ; 

LOCATE 10, 48: PRINT ",Miln": PRINT ; 

LOCA'Ili 12, 48: PRINT ",%opn VI": PRINT ; 
LOCA'Ili 13, 48: PRINT!' ",%opn V2": PRINT ; 
LOCA'Ili 15, 43: PRINT " rcs=": PRINTT ; 

LOCATE 16, 43: PRINT * yno=": PRINT ; 

LOCATE 17, 43: PRINTT "A v%=": PRINT ; 

LOCATE 18, 43: PRINT " hrd=": PRINT ; 

LOCATE 24, 50: PRINTT "P"; 

LOCATE 24, 48: PRINT CI1R$(243); 

LOCATE 24, 52: PRINT CI1R$(243); 

LOCATE 24, 43: PRINT USING MINP; 

LOCATE 24, 54: PRINT USING MAXP; 

LOCATE 25, 50: PRIN1’ "R"; 

LOCATE 25, 48: PRINT CIIR$(243); 

LOCATE 25, 52: PRINT CIIR$(243); 

LOCATE 25, 43: PRIN'T USING MINRE; 

LOCATE 25, 54: PRINT USING MAXRE; 

COLOR 14, 5 

LOCATE 15, 43: PRINT "P": PRINT ; 

LOCATE 16, 43: PRINT "R": PRINT ; 

LOCATE 17, 44: PRINT "G": PRINT ; 

LOCATE 18, 43: PRINT "C": PRINT ; 

COLOR 15, 6 

LOCATE 2, 62: PRINT "RPM/MACII LOOP": PRINT 
LOCATE 3, 64: PRINT " MANUAL ”: PRIN T ; 
LOCATE 5, 67: PRINT ",Mach": PRINT ; 

LOCATE 6, 67: PRIN T ",RPM": PRIN'T ; 

LOCATE 9, 67: PRIN T ".Mach": PRINT ; 

LOCATE 10, 67: PRIN'T ",RPM": PRIN'T ; 

LOCA'Ili 12, 67: PRINT ”,% RhsT": PRIN'T ; 
LOCATE 15, 61: PRIN'T " ach=": PRINT ; 

LOCA'TE 16, 61: PRINT " rpm=": PRINT ; 

LOCATE 24, 68: PRINT "N"; 

LOCATE 24, 66: PRINT CIIR$(243); 

LOCATE 24, 70: PRINT CIIR$(243); 

LOCATE 24, 61: PRINT USING ”####"; MINN; 
LOCATE 24, 72: PRINT USING "####"; MAXN; 
LOCATE 25, 68: PRINT "M"; 

LOCATE 25, 66: PRINT CIIR$(243); 

LOCATE 25, 70: PRINT CIIR$(243); 

LOCATE 25, 61: PRINT USING MINM; 

LOCATE 25, 72: PRINT USING "MU"; MAXM; 

COLOR 14, 6 

LOCATE 15, 61: PRINT "M": PRINT ; 

LOCATE 16, 61: PRINT "N": PRINT ; 

COLOR 15, 1 

LOCATE 5, 3: PRINT "SET POINT’: PRINT ; 
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LOCATE 9, 3: PRINT "PROCESS": PRIN T ; 
LOCATE 12, 3: PRINT "COMMAND": PRINT ; 
LOCATE 15, 3: PRINT "INPUTS": PRINT ; 
LOCATE 20, 3: PRIN T "STATUS": PRIN T ; 
LOCATE 16, 3: PRINT " elete": PRINT ; 
LOCATE 24, 3: PRINT "LIMITS"; 

COLOR 14, 1 

LOCATE 16, 3: PRIN T "D": PRIN T ; 


’Controller Program. 
1000 GOSUB 8000 


’Conversion of Digital Inputs to Engineering Units. 

PP = E(l) * XPP 
PPUSCS = PP * 14.696 
PS = E(2) * XPS 

IF PS > (PP * .999999) THEN PS = PP * .999999 
PSUSCS = PS * 14.696 

IF E(3) > 1.191 THEN 2080 ELSE 2060 

2060 'IT = 74.1826 + 105.3 * E(3) - 40.66 * E(3) A 2 + 20.54 * E(3) 3 5.21 * E(3) 4 


2080 

2090 

2100 


IT - 80.678 + 84.52 * E(3) - 12.717 * E(3) A 2 + 1.805 * E(3) A 3 - .1102 * E(3) A 4 
IE' E(4) > 1.191 THEN 2120 ELSE 2100 

TMWL = 74.183 + 105.3 * E(4) - 40.66 * E(4) A 2 + 20.54 * E(4) A 3 - 5.21 * E(4) A 4 


2120 TMWL = 80.678 + 84.52 * E(4) - 12.717 * E(4) A 2 + 1-805 * E(4) A 3 - .1102 * E(4) A 4 
2130 FRPM = E(5) * XFRPM 
PLQ = E(6) * XPLQ 
PLQUSCS = PLQ * 1 4.696 
DLP = E(7) * XDLP 


GOSUB 5000 


M = SQR(5 * (PP / PS) A .28571 - 5) 

ME' = (1 + .2 * M * M) 

KRE = 63714 * CII * M / IT A 1.4 / (ME) A 2.1 

RE = KRE * PP 

SAT = 50 + 27.34 4 PS A .296 

SAT1 = SAT * ME' 

LDPQ = (PLQ - PP) 

IF LDPQ < .5 IIIEN LDPQ = .5 
LF = CLQV * .8676 * SQR(LDPQ) 

DRPM = Ell PM 

IF DRPM < 100 THEN DRPM = 100 
KTGS = DRPM * SQR(PP) * KT / 3! / IT / LF 
TTMP = IT 

IF TTMP < 80 THEN PIMP = 80 
FKW = 100 * PP * (FRPM / 1000) A 2.26 / SQR(1TMP) 
E'B = FKW / (121 + TO / LF 
SPR = SP / 14.696 
SPR1 = SRE / KRE 

IF SPR1 > MAXSPR THEN SPR1 = MAXSPR 
IF SPR1 < MINSPR THEN SPR1 = MINSPR 
IF AUTORE = 1 THEN SPR = SPR1 
SPRU = SPR * 14.6% 

COLOR 14, 5 
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IF AUTORE = 1 TIIIZN I.OCATE 5, 43: PRINT USING "##.##"; SPRU 


’Temperature Control I,oop. 

IF 10 = l THEN 2950 
IF (TMWI. - Tl - ) > 1.5 * MXT TURN 2950 
IF FROM < 580 I'll BN 2950 
ST = ST1 

IF ST > IT 1 IIBN 2530 EI.SE 2525 
2525 IF ST < (STP - .04) TURN ST = (STP - .04) 

2530 IF AUTOT = 1 TURN 2540 ELSE 2850 

2540 IF AUTOP = 0 TURN 2610 

IF ABS(PP - SPR) < .15 TURN 2560 EI.SE 2570 
2560 IF ABS(TT - ST) < 2! TURN 2610 ELSE 2570 

2570 IF (PP / TI) > .95 * (SPR / SI) TURN 2610 ELSE 2580 

2580 IF IT > ST TURN 2610 RRSR 2590 

2590 ST = TT - .02 

2610 IF (TMWI, - ST - MX!) > 0 TURN 2620 RLSR 2650 

2620 ST = TMWI, - MXT 
FB = 0 

2650 IF ST < SAT I TURN ST = SAT1 

BT = TT - ST 

IF ABS(B'l) < .3 TURN FBI = 1 ELSE 11.1 = 0 
IF ABS(TMWL - FI) < 24 TURN FR2 = 1 EI.SE FL2 = 0 
IF BT < -5 T URN FF = 0 
IF BT > 0 TURN FT" = 1 
FBF = FB * IT * .8 

RI P = RITM1 + KIT * KTGS * DBF * BT 
IF RI P < -FBF TURN RI P = -FBF 
IF RIT > (1 - FBF) 'HI BN RIT = (1 - FBI) 

IF (TMWI, - 11) > MXT / 2 THEN 2750 El -SB 2780 
2750 LMT = 1 - (TMWI, - TP - MXT / 2) * 2 / MXT 
IF LMT < 1! / I.F TURN LMT = 1! / IT 
IF RIT > LMT TURN RPP = LMT 

2780 ALQ = KTGS * (KPT * BT + KDT * (ET - 2 * BTM1 + ETM2) / 2 / DEL) + RIT + FBF 
IF' ALQ < 0 TURN ALQ = 0 

IF ALQ > 1 TURN ALQ = 1 

IF (TMWI, - TI) > MXT / 2 TURN 2820 EI.SB 2830 
2820 IF ALQ > LMT TURN ALQ = LMT 

2830 LCMDS = ALQ * 100 
GOTO 2970 

2850 IF LCMDS > 100 THEN LCMDS = 100 

ALQ = LCMDS / 100 

IF (TMWI. - TI) > MXT / 2 TURN 2880 RUSE 2910 
2880 LMT = 1 - (TMWI, - IT - MXT / 2) * 2 / MXT 
IF LMT < 1! / I.F TURN LMT = 1! / IT 
IF AIX) > IMP TURN ALQ = IMT 
2910 RIT = ALQ 
GOTO 2970 
2950 ALQ = 0 

LCMDS = 0 

2970 IF ALQ > AI.QP + .01 THEN ALQ = ALQP + .01 

DAC(1) = ALQ 
DAC(2) = ALQ 

’Pressure - Reynolds Number Control Ixxtp. 

KPGS = 750 * KP / PP / SQR(TI) 

GF1 = 2.725 * CGV * PP / SQR('IT) 
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IF ppuscs > MAXP1 THEN PIE = 1: GOSUB 6000 
IF IE = 1 II IEN 3220 

IF ABS(RE - SRE) < .05 THEN FL7 = 1 ELSE FI-7 = 0 
H' AUTOl’ = 1 THEN 3040 El SE 3180 
FP = 1*1* - Sl’R 

IF ABS(EP) < .005 THEN 113 = 1 EI-SE M3 = 0 
RIP = RIl’Ml + El* • KI’GS * KIP » DEE 
IF RIP < 0 THEN RIP = 0 
IF RIP > 1! THEN RIP = 1! 

AG VI = KPGS * (KPP * EP + KDP * (EP - 2 * EPM1 + EPM2) 12 1 DEL) + RIP 
IF AGV1 < 0 THEN AGV1 = 0 
IF AG VI > 1! THEN AGV1 = 1! 

IF AG VI > .9 THEN AGV2 = AGV2 + .01 
IF AG VI < .7 THEN AGV2 = AGV2 - .01 
IF AGV2 < 0 THEN AGV2 = 0 
IF AGV2 > 1! THEN AGV2 = 1! 

GCMDS = AG VI * 100 
GOTO 3240 

IF GCMDS > 100 THEN GCMDS = 100 
RIP = GCMDS / 100 
AG VI = GCMDS / 100 

IF FRPM < 300 THEN LCMDS = 0 
GOTO 3240 
AG VI = 1! 

GCMDS = 100 

IF AG VI < AGV1P - .05 THEN AGV1 = AGV1P - .05 
IF AG VI > AG VI I* + .05 THEN AGV1 = AGV1P + .05 
DAC(4) = AGV1 * .8 + .2 
DAC(5) = AGV2 ♦ .8 + .2 

’Fan Speed - Mach Number Control I .nop. 

‘ IF IE = 1 THEN 3405 ELSE 3410 

3405 NCMDS1 = 0 
GOTO 3550 

3410 NCMDS1 = NCMDS 

IF FRPM > MAXN1 'll IEN NIE = 1: GOSUB 6000 
IF AUTOM = 1 THEN 3422 ELSE 3550 
3422 IF FRPM < 580 THEN AUTOM = 0: GOTO 3435 

FKM = 1 - .5 * M 

IF FKM < .71 'll IEN FKM = .71 
KMC = 520 * SQR(TT) * FKM / PP " .035 
KMGS = KMC * KMM 

IF AI.Q > .99 THEN 3432 ELSE 3440 
3432 AUTOM = 0 

NCMDS = FRPM - 500 

IF NCMDS < 0 THEN NCMDS = 0: BEEP 
3435 COLOR 14, 6 

LOCATE 3, 64: PRINT " MANUAL ": PRINT 
LOCATE 6, 62: PRINT USING "####."; NCMDS 
GOTO 3550 
3440 EM = SM - M 

IF ABS(EM) < .002 THEN FM = 1 ELSE Fl>4 = 0 
RIM = RIMM1 + KIM * EM * DEL * KMGS 

IF RIM < -100 THEN RIM = -100 

IF RIM > 100 THEN RIM = 100 
NCMDS1 = KPM * EM * KMGS + RIM + FRPM 

IF NCMDS 1 > 5450 THEN NCMDS1 = 5450 


3040 


3180 

3220 

3240 


21 


IF NCMDS1 > MAXN THEN NCMDSI = MAXN 
3550 EN = NCMDSI - FRPM 

IF ABS(EN) < 10 THEN FLN=! ELSE EEN=0 
IF ABS(EN) > 100 THEN 3560 ERSE 3570 
3560 KEN - 5 

GOTO 3575 

3570 KEN = SQR(TT) / 20 

3575 RIN = R1NM1 + EN * DEL * KIN 

IF RIN > RINMI + KEN THEN RIN = RINMI + KEN 
IF' RIN < RINMI - KEN THEN RIN = RINMI - KEN 
IF RIN < 0 THEN RIN = 0 
IF RIN > 6000 THEN RIN = 6000 
SRPM = EN * KPN + RIN 

IF SRPM > SRPM1 + 5 THEN SRPM = SRPM1 + 5 
IF SRPM < SRPM1 - 5 THEN SRPM = SRPM1 - 5 
SNRPM = SRPM / 7500 

IF SNRPM > 1! THEN SNRPM = 1! 

IF SNRPM < 0! THEN SNRPM = 0! 

DAC(6) = SNRPM 


’Fan Speed Band Warning. 

IF FRPM < 3651 THEN 3790 ERSE 3830 
3790 IF FRPM > 3549 THEN 3800 ERSE 3830 

3800 COLOR 20, 8 

LOCATE 13, 61: PRINT " SPEED BAND PRINT 
BEEP 

GOTO 3900 
3830 COLOR 14, 6 

LOCATE 13, 61: PRIN T " B : PRINT 


'Temperature Gradient Calculation. 

3900 IW = IW + 1 

IF IW a 40 THEN 3920 ELSE 3980 
3920 IF WIU > 8 THEN WLG = 0 

IF WLG < -8 THEN WLG = 0 
COLOR 14, 3 

LOCATE 20, 30: PRINT USING "###.#"; 12 * WLG 
WLG = 0 
IW = 0 

3980 WLG = WLG + TMWL - TMWI.1 


'Liquid Back Pressure Control R>op. 

SLQ = SLQSC / 14.696 

IF SLQ > 10.2 THEN SLQ = 10.2 
SLQSC = SLQ * 14.696 
ELP = PIJQ - SLQ 

IF ABS(ELP) < .4 THEN PR5 = 1 ERSE ER5 = 0 
IF ELP > .15 THEN ELP = .15 
IF ELP < -.15 'HIEN ELP = -.15 
RIL = RILM1 + ELP * KIL * DEL 
IF RIL < 0. THEN RIL = 0. 

IF RIL > 1! THEN RIL = 1! 

ALN = KPL * ELP + KDL * (ELP - ELPM1) / DEL + RIL 
IF ALN < 0. THEN ALN = 0. 

IF ALN > H THEN ALN = 1! 

DAC(3) = ALN 


22 


DAC(7) = 1 ! 

GOSUH 9000 

’Screen Pressure Drop Warning. 

13 1, PC = .47 * PP * M * M / (Ml) A 6 + .02 

IF DLP > DI.PC THEN FI/) = 1 EISE FI/> = 0 
IF M < .25 THEN H/) = 0 
IF FI/ = 1 THEN BEEP 


’Controller Screen Update. 

JD = JD + l 

IF JD = 3 THEN 4160 EUSE 4700 
4160 COLOR 14, 3 

LOCATE 21, 30: PRINT USING "###.#"; SAT1 

IF (FL1 * FI -2 * AUTOT) = 1 THEN 4190 ELSE 4200 
4190 COLOR 14, 0 

4200 LOCATE 8, 26: PRINT STRING$(6, 0) 

COLOR 14, 5 

LOCATE 21, 47: PRIN T USING ”###.##"; E12 * 20.07474 
IF F1j6 = 1 THEN COLOR 30, 5 
LOCATE 22, 48: PRINT USING "##.###"; DLP 
IF AUTORE = 0 THEN 4240 
IF (FI .7 * AUTORE) = 1 THEN 4250 FUSE 4260 
4240 IF (FI 3 * AUTOP) = 1 THEN 4250 FUSE 4260 

4250 COLOR 14, 0 

4260 LOCA'IE 8, 43: PRIN T S'I'RING$(6, 0) 

COLOR 14, 2 

LOCA1E 5, 13: PRINT USING "mi"; SLQSC 
LOCATE 9, 13: PRINT USING PLQUSCS 

LOCATE 12, 13: PRINT USING "mi"; (l - ALN) * 100 
IF FI 2) = 1 ’ll IEN 4320 EISE 4330 
4320 COLOR 14, 0 

4330 LOCATE 8, 13: PRINT STRING$(6, 0) 

COLOR 14, 6 

IF (HA * AUIOM) = 1 THEN 4360 FUSE 4370 
4360 COLOR 14, 0 

4370 LOCATE 8, 62: PRINT STRING$(6, 0) 

^(H.N-AU'TOM)^ THUN 4374 FUSE 4376 
4374 COLOR 14,0 

4376 LOCATE 11, 62: PRIN T STRINGS^, 000) 

COLOR 14, 3 

IF AUTOT = 1 THEN LOCATE 6, 26: PRINT USING ”###.#"; ST 
LOCATE 9, 26: PRINT USING "mi"; IT 
LOCATE 10, 26: PRINT USING "###.#"; TMWL 
LQ = ALQ * 100 

LOCATE 12, 26: PRINT USING "###.#"; LQ 
COLOR 14, 5 

LOCATE 9, 42: PRINT USING "###.##”; Ell * 20.07474 
GV1 = AG VI * 100 

LOCATE 12, 43: PRINT USING "###.#"; GV1 
GV2 = AGV2 * 100 

LOCATE 13, 43: PRINT USING "###.#"; GV2 
LOCATE 10, 42: PRINT USING "MUM"; RE 
COLOR 14, 6 

LOCATE 9, 62: PRINT USING "#.###"; M 
LOCATE 10, 62: PRINT USING "####."; FRPM 
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LOCATE 12, 62: PRINT USING SNRPM * 100 

JD = 0 


’Selling Previous Cycle Values. 

4700 AG VIP = AG VI 
ALQP = ALO 
EI.PM1 = LLP 
LPM2 = LPM1 
EPM1 = LP 
UTM2 - LTM1 
CTM1 = ET 
RILM1 = RIL 
RIMM1 = RIM 
RINM1 = RIN 
RIPM1 = RIP 
RITM1 = RIT 
SRPM1 = SRPM 
STP = ST 
TMWL1 = TMWL 

4900 GOTO 1000 

’Controller Input Subroutine. 

5000 AS = INKEYS 

IF AS = "" THEN 5999 
IF AS = "D" OR AS = ’’d" THEN 5950 
IF XX > 5 THEN XX = 5 
IF AS = CIIR$(13) THEN 5008 ELSE 5009 
IF XX = 5 THEN 5600 ELSE 5999 
IF J > 0 THEN 5010 ELSE 5014 

IF AS = "0" OR AS = ”1" OR AS = "2" OR A$ = "3" OR AS = "4" OR AS = "5” OR AS 
= "6" OR AS = ”7" OR AS = "8” OR AS = "9" OR AS = THEN 5100 ELSE 5999 
COLOR 16, 3 

IF AS = "T" OR AS = "t" THEN J = 1: LOCATE 15, 36: PRINT ,, ,K": PRINT ; 

IF AS = "L” OR AS = "I” THEN J = 5: LOCATE 16, 36: PRINT ",opn": PRINT ; 

COLOR 16, 5 

IF AS = "P" OR AS = ”p" THEN J = 2: LOCATE 15, 54: PRINT ",psia": PRINT ; 

IF AS = "R" OR AS = V THEN J = 4: LOCATE 16, 54: PRINT ",miln": PRINT ; 

IF AS = "G" OR AS = "g" THEN J = 6: LOCATE 17, 54: PRINT ”,opn": PRINT ; 

IF AS = "C" OR AS = "c" THEN J = 8: LOCATE 18, 54: PRINT ".in": PRINT ; 

COLOR 16, 6 

IF AS = "M" OR AS = "m" THEN J = 3: LOCATE 15, 71: PRINT ",Mach": PRINT ; 

IF AS = "N" OR AS = "n" THEN J = 7: LOCATE 16, 71: PRINT ",rpm": PRINT ; 

COLOR 16, 2 

IF AS = "B" OR AS = "b" THEN J = 9: LOCATE 15, 19: PRINT ",psi M : PRINT ; 

IF J = 1 THEN LU = 15: MU = 30 

IF J = 2 THEN LU = 15: MU = 48 

IF J = 3 THEN LU = 15: MU = 65 

IF J = 4 THEN LU = 16: MU = 48 

IF J = 5 THEN LU = 16: MU = 30 

IF J = 6 THEN LU = 17: MU = 48 

IF J = 7 THEN LU = 16: MU = 65 

IF J = 8 THEN LU = 18: MU = 48 

IF J = 9 THEN LU = 15: MU = 13 

IF J > 0 THEN BEEP 

GOTO 5999 


5008 

5009 

5010 

5014 
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51(H) 


IF J = 1 TIlliN 5105 
IF J = 2 TIIFN 5405 
IF J = 3 TIIFN 5205 
IF J = 4 TIIFN 5405 
IF J = 5 TIIFN 5405 
IF .1 = 6 TIIFN 5405 
IF J = 7 TIIFN 5350 
IF .1 = 8 TIIFN 5405 
IF J = 9 TIIFN 5105 

’Assemble Temperature and Liquid Back Pressure Inputs. 

5105 XX = XX + 1 

1)1 = 1: D2 = 2: 1)3 = 3: 1)4 = 4: 1)5 = 5 
GOSUB 5500 

AA = III * 100 + 1 12 * 10 + 113 + 115 / 10 
GO'IX) 5999 


’Assemble Mach Number Input. 


5205 


XX = XX + 1 

1)1 = 1: D2 = 3: D3 = 4: D4 = 2: 1)5 = 
GOSUB 5500 

AA = 112 / 10 + 113 / UK) + 115 / 1000 


5 


GOTO 5999 


’Assemble Fan Speed Input. 

5350 XX = XX + 1 

D1 = 1: D2 = 2: 1)3 = 3: 04 - 5: 05 = 


4 


GOSUB 5500 

AA = III * 1000 + 1 12 * 100 + 1 13 * 10 + 115 
GO'IX) 5999 


’Assemble Pressure, Reynolds Number, Injection Valve, Fxhaust Valve, and Chord Inputs. 

5405 XX = XX + 1 

01 = 1: 02 = 2: 03 = 4: 04 = 3: 05 = 5 
GOSUB 5500 

AA = III * 10 + 112 + 113 / 10 + 1 15 / 100 
GO'IX) 5999 


5500 


IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

RETURN 


J = 1 
J = 2 
J = 3 
J = 4 
J = 5 
J = 6 
J = 7 
J = 8 
J = 9 
XX = 
XX = 
XX = 
XX = 
XX = 


THEN COLOR 0, 3 


TIIFN COLOR 0, 5 
THEN COLOR 0, 6 
THEN COLOR 0, 5 
THEN COLOR 0, 3 
THEN COLOR 0, 5 
THEN COLOR 0, 6 
THEN COLOR 0, 5 

Dl^IENHf^VAUAS): LOCATE LU, MU + Dl: PRINT USING HI 

D2 H EN 2 = VAIUAS : LOCATE LU, MU + D2: PRINT USING H2 

D3 THEN 113 = VAIiA$): LOCATE LU, MU + 1)3: PRINT USING H3 

D4 THEN LOCATE LU, MU + 1)4: PRINT CHR$(46); 

D5 THEN 115 = VAL(A$): LOC3ATE LU, MU + D5: PRINT USING "# ; 115 
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’Set Temperature Input. 

5600 IF J = 1 TI IEN 5620 ELSE 5660 

5620 STI = AA 

IF STI > MAXT THEN STI = MAXI 
IF STI < MINT THEN STI = MIN I’ 
IF STI < SAT! THEN STI = SAT1 
COLOR 14, 3 

LOCATE 3, 26: PRINT " AUTO PRIN T ; 


AUTOT = 1 

LOCATE 5, 26: PRTNT USING "###.#"; STI 


’Set Pressure Input. 

5660 IF J = 2 THEN 5670 ELSE 5730 

5670 SP = AA 

IF SP < MINI 1 THEN SP = MINP 
IF SP > MAXP THEN SP = MAXP 
COLOR 14, 5 

LOCATE 3, 46: PRINT " AUTOP PRINT ; 
LOCATE 6, 43: PRINT " ": PRINT ; 

LOCATE 5, 43: PRINT USING "##.##"; SP 
AUTOP = 1 
AUTORE = 0 

’Set Mach Number Input. 

5730 IF J = 3 THEN 5740 ELSE 5750 

5740 SM = AA 

IF’ SM < MINM THEN SM = MINM 
IF SM > MAXM THEN SM = MAXM 
COLOR 14, 6 

LOCATE 3, 64: PRINT " AUTO ": PRINT ; 

LOCATE 6, 62: PRIISTI 

AUTOM = 1 

LOCATE 5, 62: PRINT USING "#.###"; SM 


’Set Reynolds Number Input. 

5750 IF J = 4 THEN 5760 ELSE 5770 

5760 SRE = AA 

IF SRE > MAXRE THEN SRE = MAXRE 
IF SRE < MINRE THEN SRE = M1NRE 
COLOR 14, 5 
AUTORE = 1 
AUTOP = 1 

LOCATE 6, 43: PRINT USING "##.##"; SRE 
LOCATE 3, 46: PRINT " AUTORE *: PRIN T ; 


’Set Liquid Injection Valve Input. 

5770 IF J = 5 THEN 5780 ELSE 5810 

5780 LCMDS = AA 
COLOR 14, 3 

LOCATE 3, 29: PRINT " MANUAL PRINT ; 
LOCATE 5, 26: PRINT " 

LOCATE 6, 26: PRINT " 

AUTOT - 0 


’Set Gas Exhaust Valve Input. 

5810 IF J = 6 THEN 5820 EI^E 5860 

5820 GCMDS = AA 
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COLOR 14, 5 
LOCATE 3, 46: PRINT 
AUTOP = 0 
AUTORE = 0 
LOCATE 5, 43: PRINT 
I .OCA'IE 6, 43: PRINT 


MANUAL ": PRINT ; 


■: PRINT ; 
": PRINT ; 


’Set Pan Speed Inpul. 

5860 IF J = 7 THEN 5870 ELSE 5935 

5870 NCMDS = AA 

IF NCMDS > MAXN THEN NCMDS = MAXN 
COLOR 14, 6 
LOCATE 5, 62: PRINT ■ 

LOCATE 3, 64: PRINT " MANUAL PRINT ; 

LOCATE 6, 62: PRINT USING "####."; NCMDS 
AUTOM = 0 

’Set Chord Ixngth Input. 

5935 IF J = 8 THEN 5937 ELSE 5943 

5937 CHIN = AA 

IF CHIN < M1NCIIIN THEN CHIN = MINC1I1N 
IF CHIN > MAXCIIIN THEN CIIIN = MAXCHIN 
COLOR 14, 5 

LOCATE 20, 48: PRIN T USING "##.##”; CIIIN 
CII = CIIIN * .0254 

’Set Liquid Back Pressure Inpul. 

5943 IF J = 9 THEN 5944 ELSE 5950 

5944 SLQSC = AA 

IF SLQSC > MAXLQSC 'IlIEN SLQSC = MAXLQSC 
COLOR 14, 2 

LOCATE 5, 13: PRINT USING "###.#”; SLQSC 

5950 II- .1 = 0 THEN 5962 

IF J = I IlIEN COLOR 14, 3 
IF J = 2 THEN COLOR 14, 5 
IE J = 3 'IlIEN COLOR 14, 6 
IF J = 4 THEN COLOR 14, 5 
IF J = 5 THEN COLOR 14, 3 
IF J = 6 THEN COLOR 14, 5 
IF J = 7 THEN COLOR 14, 6 
IF J = 8 THEN COLOR 14, 5 
IF J = 9 THEN COLOR 14, 2 
LOCATE LU, MU + 1: PRINT " ": PRIN T ; 


5962 IF J > 0 THEN BEEP 

IF IE = 1 THEN COLOR 0, 6 

IF IE = 1 THEN BEEP: LOCATE 20, 60: PRINT " PRINT 

IF PIE = 1 THEN LOCATE 21, 60: PRINT " PIE = 0 

IF NIE = 1 THEN LOCATE 21, 60: PRINT " NIE = 0 

J = 0 
XX = 0 
IE = 0 

5999 RETURN 
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Emergency Stop Subroutine. 

6000 IB = 1 

COLOR 20, 8 

LOCATE 20, 60: PRINT " EMERGENCY STOP PRINT ; 

COLOR 4, 8 

IF PIE = I THEN LOCATE 21, 60: PRINT " Pressure Limit M 

IF NIB = 1 THEN LOCATE 21, 60: PRINT! 1 M Fan RPM Limit M 

COLOR 4, 8 
AUTOT = 0 
AUTOP = 0 
AUTOM = 0 
AUTORE = 0 

LOCATE 3, 29: PRINT H MANUAL PRINT ; 

LOCATE 3, 46: PRINT " MANUAL PRINT ; 

LOCATE 3, 64: PRIN T " MANUAL PRINT ; 

RETURN 

’Analog to Digital Conversion Input Subroutine. 

8000 OUT COMMAND.REGISTER, ZSTOP 
ZEMP = INP(ZATA.REGISTER) 

WAI T ZTATUS.REGISTER, ZOMMAND.WAIT 
OUT ZOMMAND.REGISTER, ZCLEAR 
ZDGAIN = 1 
FOR I = 1 TO 7 
ZDCIINL =1-1 

WAIT ZTATUS.REGTS TER, WRITE. WAIT, WRTTE.WAIT 

WAIT ZTATUS. REGISTER, ZOMMAND.WAI T 

OUT ZOMMAND.REGISTER, ZADIN 

WAIT ZTATUS.REGISTER, WRTTE.WAIT, WRTTE.WAIT 

OUT ZATA.REGISTER, ZDGAIN 

WAIT ZTATUS.REGISTER, WRITE. WAIT, WRTTE.WAIT 

OUT ZATA.REGISTER, ZDCIINL 

WATT ZTATUS.REGISTER, READ.WAIT 

ZOW = INP(ZATA.REGISTER) 

WATT ZTATUS.REGISTER, READ.WAIT 
ZIGII = INP(ZATA. REG ISTER) 

ZOLT# = ZIGII * 256 + ZOW 

IF ZOLT# > 32767 17IEN ZOLT # = ZOLT# - 65536! 

WATT ZTATUS.REGISTER, ZOMMAND.WAIT 
ZTATUS = INP(ZTATUS. REGISTER) 

IF (ZTATUS AND &II80) THEN GOTO 8000 
ZOL# = 5 * ZOLT# / 32768 
E(I) = ZOL# 

NEXT I 


’Sensor Failure Detection. 

El 1 = 0(1): IF El 1 <0 THEN El 1 = 0 

IF E(l) < .43 THEN E(l) = .43: FL10 = 1 
E12 = E(2): IF E12 < 0 THEN E12 = 0 

IF E(2) < .4 THEN E(2) = .4: FL10 = 2 
IF E(3) < 0 THEN E(3) = 0!: FL10 = 3 
IF E(4) < 0! THEN E(4) = 0: FL10 = 4 
IF E(5) < .001 IIIEN E(5) = .001 
IF E(6) < .3 IIIEN E(6) = .3 
IF ABS(E(7)) > I ! THEN FL10 = 7 
IF I *1.10 >= 1 THEN 8340 EI-SE 8380 
8340 COLOR 4, 8 
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LOCATE 21, 60: PRINT " Sensor Failure PRINT ; 

IF FI. 10 = 1 THEN LOCATE 22, 60: PRINT " Total Pressure " 

IF FI, 10 = 2 THEN LOCATE 22, 60: PRINT " Static Pressure " 

IF IT10 = 3 THEN LOCATE 22, 60: PRINT " Gas Temperature " 

IF FL10 = 4 THEN LOCATE 22, 60: PRINT " Wall Temperature" 

IF FL10 = 7 THEN LOCA TE 22, 60: PRIN T " Screen Pressure " 

GOSUB 6000 
GOTO 8400 
8380 COLOR 6, 6 

LOCATE 21, 60: PRINT " PRINT ; 

LOCATE 22, 60: PRINT " PRINT ; 

8400 FL10 = 0 
RETURN 


’Digital 

to Analog Conversion Output Subroutine. 

9000 

FOR IK = 1 TO 7 

DR1 = INT(DAC(IK) * 4095) 

III = INT(DRI / 16) 

DR2 = INT(DR1 - III * 16) 

Il< IK < 3 THEN 9050 FUSE 9070 

9050 

LO = INT(DR2 * 16) + 2 * IK - 2 
GO TO 9080 

9070 

LO = INT(DR2 * 16) + 2 *IK - 1 

9080 

BEMP = INP(&I 1225) 

IF HEMP = 0 THEN 9110 

9100 

GOTO 9080 

9110 

OUT &II224, LO 

9120 

CEMP = INP(&H225) 

IF CEMP = 16 THEN 9150 
GOTO 9120 

9150 

OUT & 1 1224, III 
NEXT IK 
RETURN 


10000 


END 


Appendix B 

Controller Simulator Source Code Listing 


’0.3-m Transonic Cryogenic Tunnel Controller Simulator. 
CIJS 
Cl ALAR 


' Program Const ants. 

DEL = .1 
CGV = 8 
C1T = .18 

CHIN = .18 * .0245 
CLQV = 4 
IW = 39 
K = 1 
KDL = 0! 

KDP = 0! 

KDT = 0! 

KIL = .02 
KIM = 0! 

KIN = .4 
KIP = .05 
KIT = .1 
KMM = .3 
KP = 3 
KPL * .2 
KPM = 4.5 
KPN = .6 
KPP = 1 ! 

KPT = 1 ! 

KT = .04 
MXT = 40 
PIjOI = 3! 

PLQ = 3! 

PP = 1 ! 

PP1 = 1! 

SLQSC = 17! 

STP = 300! 

TMWL = 300! 
TMWL1 = 300! 

IT = 300! 

TT1 = 300! 


’Maximum - Minimum Limits Screen layout. 

COLOR 14, 9 
CLS 

LOCATE 3, 16: PRINT " Maximum - Minimum Limits for 0.3-m TCT Controller " 
COLOR 12, 9 

LOCATE 5, 25: PRINT ’Tunnel” 

LOCATC 6, 24: PRINT "Maximun" 

LOCATC 5, 38: PRINT "Set" 

LOCATE 6, 36: PRINT "Maximun" 

LOCATC 5, 60: PRINT "Set" 

LOCATC 6, 58: PRINT "Minimum" 

LOCATC 5, 72: PRINT "Tunnel" 
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LOCATE 6, 71: PRINT "Minimum 


COLOR 7, 9 

LOCATE 7, 3: PRINT ’Temperature" 

LOCA1T 8, 4: PRINT "(kclvin)" 

LOCATE 10, 3: PRIN T "Pressure" 

LOCATE 11, 4: PRINT "(psia)" 

LOCATE 13, 3: PRIN T "Normal I’an Speed” 

LOCATE 14, 4: PRIN T "(rpm)" 

LOCATE 16, 3: PRINT "Mach Number" 

LOCATE 19, 3: PRINT "Reynolds Number" 

LOCATE 20, 4: PRINT "(Miln)" 

LOCATE 22, 3: PRINT "Accept Tunnel Max./Min. Limits and Start 0.3-m TCT Controller. 


LOCATE 


23, 3: PR [N'T "Accept Set Max./Min. Limits and Start 0.3-m TCT Controller." 

COLOR 15, 9 

LOCATE 7, 51: PRIN T T’ 

LOCATE 7, 49: PRINT CIIR$(243) 

LOCATE 7, 53: PRIN T CI1R$(243) 

LOCATE 10, 51: PRIN T "P" 

LOCATE 10, 49: PRINT CI1R$(243) 

LOCATE 10, 53: PRIN T CI1R$(243) 

LOCATE 13, 51: PRIN T "N" 

LOCATE 13, 49: PRINT CIIR$(243) 

LOCATE 13, 53: PRINT CIIR$(243) 

LOCATE 16, 51: PRIN T ”M” 

LOCATE 16, 49: PRINT CMR$(243) 

LOCATE 16, 53: PRIN T CI!R$(243) 

LOCATE 19, 51: PRINT ”R" 

LOCATE 19, 49: PRINT CIIR$(243) 

LOCATE 19, 53: PRINT CIIR$(243) 


COLOR 14, 9 
LOCATE 7, 3: PRINT "T" 
LOCATE 10, 3: PRINT "P" 
LOCATE 13, 3: PRINT "N" 
LOCATE 16, 3: PRIN T "M" 
LOCATE 19, 3: PRIN T "R" 
LOCATE 22, 3: PRINT "A" 
LOCATE 23, 10: PRIN T ”S” 


Tunnel Maximum - Minimum Limits. 
MAXT1 = 340! 

MINT1 = 77! 

MAX PI = 88! 

MINPl = 14.696 
MAXN1 = 5600! 

MINN1 = 0! 

MAXM1 = .995 
M1NM1 = .15 
MAXRE1 = 50! 

MINREl = 1! 

MAXCIIIN = 15.75 
MINCIIIN = .394 
MAXCII = .4 
MINCII = .01 
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MAXI.QSC = 150! 

MAXSI'R = MAXIM / 14.6% 

MINSPR = MINP1 / 14.696 

MAXT = MAXT1 
MINT = MINT1 
MAXP = MAXP1 
MINP = MINP1 
MAXN = MAXN1 
MINN = MINN l 
MAXM = MAXM1 
MINM = MINM1 
MAX RE = MAXRE1 
MINRE = MINRE1 

COLOR 10, 9 

LOCATE 7, 25: PRINT USING "###.#"; MAXTl 
LOCATE 7, 37: PRINT USING "###.#”; MAXTl 
LOCATE 10, 25: PRINT USING "##.##"; MAXIM 
LOCATE 10, 59: PRINT USING "##.##”; MINP1 
LOCATE 13, 25: PRINT USING "####."; MAXN1 
LOCATE 13, 60: PRINT USING "####.”; MINN1 
LOCATE 16, 25: PRIN T USING ”#.###"; MAXM1 
LOCATE 16, 60: PRINT USING "#.###"; MINM1 
LOCATE 19, 25: PRINT USING "##.##"; MAXRE1 
LOCATE 19, 60: PRINT USING "##.##"; MINRE1 
LOCATE 7, 69: PRINT "Saturation" 

LOCATE 8, 69: PRINT "Temperature" 

LOCATE 10, 71: PRINT "Ambient" 

LOCATE 11, 71: PRINT "Pressure" 

LOCATE 13, 73: PRIN T USING "####."; MINN1 
LOCATE 16, 73: PRINT USING "#.###"; M1NM1 
LOCATE 19, 73: PRINT USING "#0#"; MINREI 


’Maximum - Minimum Limits Inputs. 

5 AS = INKEY$ 

IF AS = "T" OR AS = "1" THEN GOSUB 10 
IF AS = "P" OR AS = "p" THEN GOSUI3 100 
IF AS = "N" OR AS = "n" THEN GOSUB 2(X) 
IF AS = "M” OR AS = "m" THEN GOSUB 300 
IF AS = "R” OR AS = "r" THEN GOSUB 400 
IF AS = "A" OR AS = "a" THEN 500 
IF AS = "S" OR AS = "s" THEN 600 

GO TO 5 


’'Temperature Minimum Limit Subroutine. 

10 COLOR 15, 9 

BEEP 
LR = 7 
LC = 58 

LOCATE LR, LC + 1: PRINT " ,K" 

XX = 0 

20 B$ = INKEYS 

IF B$ = "D" OR BS = "d" '11 1EN 30 ELSE 40 

30 LOCATE LR, LC + 1: PRINI 

RETURN 

40 IF BS = "0" OR BS = ”1” OR BS = ”2" OR BS = ”3" OR BS = "4” OR BS = "5” OR BS 
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"6" OR B$ = "7" OR 13$ = "8" OR 13$ = "9" OR 13$ = n ." THEN 50 EESE 20 
SO XX = XX + 1 

II' XX = 1 THEN 111 = VAISES): EOCATE ER, EC + 1: PRINT USING HI 

IP' XX = 2 THEN 112 = VAI<H$): EOCATE ER, EC + 2: PRINT USING 112 

IP' XX = 3 THEN 1 13 = VAUBS): EOCATE ER, I.C + 3: PRINT USING H3 

IE XX = 4 THEN EOCATE ER, EC + 4: PRINT CIIR$(46) 

IE XX = 5 THEN 115 = VAI.(B$): EOCATE ER, EC + 5: PRINT USING 115 

IP' XX = 5 THEN 60 EESE 20 
60 AA = III * 100 + 112 * 10 + 113 + 115 / 10 
MINI’ = AA 

IE MINT >= MAXT1 IIIEN MIN’!’ = MAXT1 
IE MINT <= MINT1 THEN MINI’ = MINT1 
70 CS = INKEY$ 

IF CS = CIIR$(I3) THEN 80 ELSE 70 

80 BEEP 

COEOR 14, 9 

EOCAIE ER, EC + 1: PRINT USING "###.# MINT 
RETURN 


’Pressure Maximum Eimit Subroulinc. 

100 COEOR 15, 9 
BEEP 
EC = 36 
ER = 10 

EOCATE ER, EC + 1: PRINT " ,psia" 

XX = 0 

110 B$ = INKEYS 

IP' B$ = "D" OR B$ = "d" THEN 120 EI-SE 130 
120 EOCATE ER, EC + 1: PR INI 

130 RlirUR N , = „ () „ QK H$ = „ r QR H$ _ "2" OR B$ = "3" OR B$ = "4” OR 13$ = "5" OR 13$ 

"6" OR 13$ = "7” OR 13$ = "8" OR 13$ = ”9" OR 13$ = THEN 140 ElisE 110 
140 XX = XX + 1 

IE XX = 1 'IIIEN 111 = VAI>(B$): EOCATE ER, EC + 1: PRINT USING III 

IE XX = 2 THEN 112 = VAI^BS): EOCATE ER, EC + 2: PRINT USING H2 

IE XX = 4 'HIEN 113 = VAI<B$): EOCA'I’E ER, LC + 4: PRINT USING ID 

IE XX = 3 THEN EOCAIE ER, EC + 3: PRINT CIIR$(46) 

IF XX = 5 THEN 115 = VAI<B$): EOCATE ER, EC + 5: PRINT USING H5 

IF XX = 5 THEN 150 EESE 110 
150 AA = III * 10 + 112 + 113 / 10 + 115 / 100 
MAXP = AA 

IF MAXP <= MINIM THEN 160 EI^E 170 

160 MAXP = MINIM ^ w . 

COEOR 20, 8: EOCAIE ER + 1, 36: PRINT " Warning! Set Max. = Tunnel Min. 

Hi) = 1 
GOTO 180 

170 COEOR 14, 9: EOCAIE ER + 1, 36: PRINT " 

EIjO = 0 

180 IF MAXP > MAXIM IIIEN MAXP = MAXP1 

190 C$ = INKEYS 

IE CS = CHR$(13) THEN 195 EHSE 190 

195 BEEP 

COEOR 14, 9 

EOCATE ER, EC + 1: PRINT USING "##J# MAXP 
RETURN 
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'Normal Fan Spew) Maximum r.imil Subroutine. 

2(H) COLOR 15, 9 

beep 

I.C = 36 
LR = 13 

LOCATE l.R, LC + 1: PRINT " ,rpm" 

XX = 0 

210 B$ = INKEYS 

II-' BS = "O'* OR 11$ = M" THEN 220 EliJE 230 
220 LOCATE LR, LC + 1: PRINT " " 

RETURN 

23() IF B$ = "0" OR B$ = ”\" OR B$ = "2" OR B$ = "3" OR B$ = "4" OR B$ = "5" OR BS 

"6" OR B$ = "7” OR BS = "8" OR BS = "9" OR BS = V THEN 240 EI.SE 210 
240 XX = XX + 1 

IF XX = 1 THEN HI = VAI^BS): LOCATE LR, I.C + 1: PRINT USING HI 

IF' XX = 2 r HIEN 112 = VAI<B$): LOCATE LR, LC + 2: PRINT USING H2 

IF XX = 3 THEN 113 = VAI<BS): LOCATE LR, LC + 3: PRINT USING H3 

IF XX = 5 THEN LOCATE LR, LC + 5: PRINT CIIRS(46) 

IF XX = 4 THEN 115 = VAISES): LOCATE LR, LC + 4: PRINT USING 115 

IF' XX = 5 IHEN 250 ELSE 210 
250 AA = IT1 * 1000 + 112 * 100 + 113 * 10 + 115 
MAXN = AA 

IF' MAXN <= M1NN1 THEN 260 ELSE 270 
260 MAXN = MINN1 

COLOR 20, 8: LOCATE LR + 1, 36: PRINT " Warning! Set Max. = Tunnel Min. " 

I-TjO = 1 
GOTO 280 

270 COLOR 14, 9: LOCATE LR + 1, 36: PRINT " 

FIjO = 0 

280 IF MAXN > MAXN1 THEN MAXN = MAXN1 

290 C$ = INKEYS 

IF C$ = C1IR$(13) THEN 295 ELSE 290 

295 BEEP 

COLOR 14, 9 

LOCATE LR, LC + 1: PRINT USING "####. MAXN 
RETURN 


’Mach Number Maximum Limit Subroutine. 

300 COLOR 15, 9 
BEEP 
LC = 36 
LR = 16 

LOCATE LR, LC + 1: PRINT " ,Mach” 

XX = 0 

310 BS = INKEYS 

IF BS = "D" OR BS = "d" THEN 320 ELSE 330 

320 LOCATE LR, IT: + I: PRIN1 

RETURN 

330 IF’ BS = "0" OR BS = "1" OR BS = "2" OR BS = "3" OR B$ = H 4" OR BS = "5" OR BS 

"6" OR BS = "7" OR BS = "8" OR BS = "9" OR BS = THEN 340 ELSE 310 
340 XX = XX + 1 

IF XX = 1 THEN III = VAI^BS): LOCATE LR, LC + 1: PRINT USING III 

IF' XX = 3 THEN 112 = VAI^BS): LOCATE LR, LC + 3: PRINT USING H2 

IF XX = 4 'IHEN 1 13 = VAI^BS): LOCATE: LR, LC + 4: PRINT USING 113 

IF XX = 2 THEN LOCATE I.R, LC + 2: PRINT CHR$(46) 

IF XX = 5 IHEN 1 15 = VAL(B$): LOCATE LR, LC + 5: PRINT USING FIS 

IF XX = 5 IHEN 350 ELSE 310 
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350 AA = III + 112 / 10 + 113 / 100 + 115 / 1000 
MAXM = AA 

IF MAXM <= MINM1 THEN 360 EI5E 370 

360 MAXM = MINM1 w ^ _ „ 

COLOR 20, 8: LOCATE LR + 1, 36: PRINT " Warning! Set Max. = Tunnel Min. 

FID = I 
GOTO 380 

370 COLOR 14, 9: LOCATE LR + 1, 36: PRINT " 

FID = 0 

380 IF MAXM > MAXM1 THEN MAXM = MAXM1 

390 C$ = INKEYS 

IF C$ = CIIR$(13) THEN 395 EI5E 390 

395 BEEP 

COLOR 14, 9 

LOCATO LR, LC + 1: PRINT USING "#.### MAXM 
RETORN 

’Reynolds Number Maximum Limit Subroutine. 

400 COLOR 15, 9 
BEEP 
LC = 36 
LR = 19 

LOCATE LR, LC + 1: PRINT " ,Miln” 

XX = 0 

410 B$ = INKEYS 

IF US = "D" OR US = "d" WEN 420 EI5E 430 
420 LOCATE LR, LC + l: PRIN T " " 

430 IF B$ = "0" OR BS = "1” OR B$ = "2" OR US = "3" OR US = "4" OR US = "5" OR BS 

"6" OR BS = ”7" OR BS = ”8" OR US = "9” OR US = V THEN 440 EI5E 410 

440 X X X 4* j 

IF XX = 1 'THEN III = VAL(B$): LOCATE LR, LC + 1: PRINT USING III 

IF XX = 2 THEN 112 = VAIiBS): LOCATE LR, LC + 2: PRINT USING H2 

IF XX = 4 THEN 113 = VAI^BS): LOCATC LR, LC + 4: PRINT USING ID 

IF XX = 3 THEN LOCATE LR, LC + 3: PRINT CIIR$(46) 

IF XX = 5 TOEN 115 = VAI^BS): LOCATE I.R, I.C + 5: PRIN T USING H5 

IF XX = 5 THEN 450 EI5E 410 
450 AA = HI * 10 + 112 + 113 / 10 + 115 / 100 
MAXRE = AA 

IF MAXRE <= MINRE1 THEN 460 EI5E 470 
460 MAXRE = MINRE1 _ „ 

COLOR 20, 8: LOCATE LR + 1, 36: PRIN T " Warning! Set Max. = Tunnel Min. 

FID = 1 

GOTO 480 

470 COLOR 14, 9: LOCATE LR + 1, 36: PRIN T " 

FID = 0 

480 IF MAXRE > MAXRE1 TOEN MAXRE = MAXRE1 

490 C$ = INKEYS 

IF C$ = C1IR$(13) TOEN 495 EI5E 490 

495 BEEP 

COLOR 14, 9 

LOCATE LR, LC + 1: PRIN T USING ”##.## MAXRE 
RETORN 

’Starts the Controller with the Tunnel Maximum - Minimum Limits. 

500 XX = 0 

GOTO 800 
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'Starts the Controller with the Set Maximum - Minimum Limits. 
600 IF Mi) = I THEN GOTO 5 

XX = 0 

MAXSPR = MAXP / 14.696 
MINSPR = MINP / 14.6% 

GOTO 800 


’Controller Simulation. 
’Controller Screen layout. 
800 COLOR 0, 0 

cfs 


COLOR 15, 4 

LOCATE 1, 26: PRINT " 0.3-m 'ITJNNEL T-P/R-M CONTROLLER " 

FOR I = 2 TO 25 
COLOR 15, 1 

LOCATE I, 2: PRINT STRING$(10, 0); 

IF’ I = 4 THEN LOCATE 4, 2: PRINT STRING$(10, 205); 

IF I = 14 THEN LOCATE 14, 2: PRINT S'1R)NG$(10, 205); 
IF I = 19 THEN LOCATE 19, 2: PRINT STRING$(10, 205); 
IF I = 23 THEN LOCATE 23, 2: PRINT STRING$(18, 205); 
COLOR 15, 2 

LOCATE I, 12: PRINT STRING$(12, 0); 

IF I = 4 THEN LOCATE 4, 12: PRINT STR1NG$(12, 205); 
IF I = 14 THEN LOCATE 14, 12: PRINT STRING$(12, 205): 
IF I = 19 THEN LOCATE 19, 12: PRINT STRING$(12, 205): 
IF I = 23 THEN LOCATE 23, 12: PRINT STRING$(18, 205); 
COLOR 15, 3 

LOCATE I, 24: PRINT STRING$(I7, 0); 

IF I = 4 THEN LOCATE 4, 24: PRINT STRING$(17, 205); 

IF' I = 14 THEN LOCATE 14, 24: PRINT STRING$(17, 205) 
IF I = 19 THEN LOCATE 19, 24: PRINT STRING$(17, 205) 
IF I = 23 THEN LOCATE 23, 24: PRINT STRING$(18, 205) 
COLOR 15, 5 

LOCATE I, 41: PRINT STRING$(18, 0); 

IF I = 4 THEN LOCATE 4, 41: PRINT STRING$(18, 205); 

IF I = 14 THEN LOCATE 14, 41: PRINT STRING$(18, 205) 
II- I = 19 THEN LOCATE 19, 41: PRINT STRING$(18, 205) 
IF I = 23 THEN LOCATE 23, 41: PRINT STOING$(18, 205) 
COLOR 15, 6 

LOCATE 1, 59: PRINT STRING$(18, 0); 

IF I = 4 THEN LOCATE 4, 59: PRINT STRING$(18, 205); 

IF I = 14 THEN LOCATE 14, 59: PRINT STRING$(t8, 205) 
IF I = 19 THEN LOCATE 19, 59: PRINT STRING$(18, 205) 
IF I = 23 THEN LOCATE 23, 59: PRINT STRING$(18, 205) 

NEXT I 


COLOR 15, 2 

LOCATE 2, 14: PRINT " LN PUMP PRINT ; 
LOCATE 3, 14: PRINT ” AUTO ”: PRINT ; 
LOCATE 5, 18: PRINT ”,psia": PRINT ; 
LOCATE 9, 18: PRINT ",psia": PRINT ; 
LOCATE 12, 18: PRINT ",%opn": PRINT ; 
LOCATE 15, 12: PRINT " =": PRINT ; 
LOCATE 5, 13: PRINT USING "###.#"; SI.QSC 
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COLOR 14, 2 

I .OCA'Ili 15, 12: PRINT "11": PRINT ; 

COLOR 15, 3 

LOCATE 2, 28: PRINT TUMP LOOP": PRINT ; 
LOCATE 3, 29: PRINT " MANUAL PRINT ; 
LOCATE 5, 31: PRIN T \K(T'inal)": PRINT ; 

LOCATE 6, 31: PRINT ",K(Usc)": PRINT ; 

LOCATE 9, 31: PRINT ",K-GN2": PRINT ; 

LOCATE 10, 31: PRINT ",K-WALL": PRINT ; 

LOCATE 12, 31: PRINT ",%opn": PRINT ; 

LOCATE 15, 25: PRIN T " cmp=": PRINT ; 

LOCATE 16, 25: PRINT "A Q%=": PRINT ; 

LOCATE 20, 25: PRINT "GRAD= ,K/mt": PRINT ; 
LOCATE 21, 25: PRINT " SAT= ,K'': PRINT ; 
LOCATE 24, 32: PRINT "T”; 

LOCATE 24, 30: PRINT CIIR$(243); 

LOCATE 24, 34: PRIN T CIIR$(243); 

LOCATE 24, 26: PRINT USING ”###"; MI N'T; 

LOCATE 24, 36: PRINT USING "###"; MAXI’; 

COLOR 14, 3 

LOCATE 15, 25: PRIN T "T": PRINT ; 

LOCATE 16, 26: PRINT "L": PRINT ; 

COLOR 15, 5 

LOCATE 2, 45: PRINT "Pl/Rc LOOP": PRINT ; 
LOCATE 3, 46: PRINT " MANUAL ": PRINT ; 
LOCATE 5, 48: PRIN T ",Psia": PRINT ; 

LOCATE 6, 48: PRIN T ",Miln": PRINT ; 

LOCATE 20, 42: PRIN T "CIIORD= 7.09, in": PRINT ; 
LOCATE 9, 48: PRIN T H ,Psin H : PRINT ; 

LOCATE 21, 42: PRIN T "P st= ,psia": PRINT ; 

LOCATE 22, 42: PRINT "Del P= ,psi": PRINT ; 

LOCATE 10, 48: PRINT ”,Miln": PRINT ; 

LOCATE 12, 48: PRINT ",%opn VI": PRINT ; 
LOCATE 13, 48: PRINT ",%opn V2": PRINT ; 
LOCATE 15, 43: PRINT " rcs=": PR I N'T ; 

LOCATE 16, 43: PRIN T " yno=": PRINT ; 

LOCATE 17, 43: PRINT "A v%=": PRINT ; 

LOCATE 18, 43: PRIN T " hrd=”: PRINT ; 

LOCATE 24, 50: PRINT "P"; 

LOCATE 24, 48: PRIN T CIIR$(243); 

LOCATE 24, 52: PRIN T CIIR$(243); 

LOCATE 24, 43: PRIN T USING MINP; 

LOCATE 24, 54: PRINT USING MAXP; 

LOCATE 25, 50: PRINT "R"; 

LOCATE 25, 48: PRINT CIIR$(243); 

LOCATE 25, 52: PRINT CIIR$(243); 

LOCATE 25, 43: PRINT USING M1NRE; 

LOCATE 25, 54: PRINT USING MAXRE; 

COT. OR 14, 5 

LOCATE 15, 43: PRINT "P": PRINT ; 

LOCATE 16, 43: PRINT "R": PRINT ; 

LOCATE 17, 44: PRINT "G": PRINT ; 
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LOCATE 18, 43: PRINT "C": PRINT ; 


COLOR 15, 6 

LOCATE 2, 62: PRINT "RPM/MACII I, OOP": PRINT 
LOCATE 3, 64: PRINT " MANUAL PRINT ; 
LOCATE 5, 67: PR INI’ ",Mach": PRINT ; 

LOCATE 6, 67: PRINT ",RPM": PRINT ; 

LOCATE 9, 67: PRINT ",Mach": PRINT ; 

LOCATE 10, 67: PRINT ",RPM": PRIN T ; 

LOCATE 12, 67: PRINT ",% RhsT": PRINT ; 
LOCATE 15, 61: PRINT " ach=": PRINT ; 

LOCATE 16, 61: PRINT ■ rpm=": PRINT ; 

LOCATE 24, 68: PRINT "N"; 

LOCATE 24, 66: PRINT CIIR$(243); 

LOCATE 24, 70: PRINT CIIRS(243); 

LOCATE 24, 61: PRINT USING "####"; MINN; 
LOCATE 24, 72: PRINT USING "####"; MAXN; 
LOCATE 25, 68: PRINT "M"; 

LOCATE 25, 66: PRINT CIIR$(243); 

LOCATE 25, 70: PRINT CIIR$(243); 

LOCATE 25, 61: PRINT USING MINM; 

LOCATE 25, 72: PRINT USING ".###"; MAXM; 

COLOR 14, 6 

LOCATE 15, 61: PRINT "M": PRINT ; 

LOCATE 16, 61 : PRINT "N": PRINT ; 

COLOR 15, 1 

LOCATE 5, 3: PRINT "SET POINT": PRIN T ; 
LOCATE 9, 3: PRIN T "PROCESS": PRINT ; 

LOCATE 12, 3: PRINT "COMMAND": PRINT ; 
LOCATE 15, 3: PRINT "INPU’l'S": PRINT ; 

LOCATE 20, 3: PRINT "STATUS": PRINT ; 

LOCATE 16, 3: PRIN T " eletc": PRINT ; 

LOCATE 24, 3: PRINT "LIMITS"; 

COLOR 14, I 

LOCATE 16, 3: PRINT "D": PRINT ; 


’Simulalion Controller Program. 

1000 CRPM = FRPM 

IE FRPM < 0 THEN FRPM = 0 
IE CRPM < 50 THEN CRPM = 50 
CRM = 2.77 - CRPM / 157 / SQR(1T) * PP A .035 
IF CRM < .001 THEN CRM = .001 
M = 1.67 - SQR(CRM) 

IF M > .999 'THEN M = .999 
PS = PP / (1 + .2 * M A 2) A 3.5 
PPUSCS = PP * 14.696 
PSUSCS = PS * 14.696 
PLQUSCS = PLQ * 14.6% 


GOSUB 5000 

MF = (1 + .2 * M * M) 

KRE = 63714 * CII * M / IT A 1.4 / (MF) " 2.1 
RE = KRE * PP 
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SAT = 50 + 27.34 * PS * .2% 

SAT1 = SAT * MI 7 
LDPQ = (PTO - PP) 

IF LDPQ < .5 THEN LDPQ = .5 
LF = CI.QV * .8676 * SQR(LDPQ) 

DRPM = FRPM 

IF DRPM < 100 THUN DRPM = 100 
KTGS = DRPM * SQR(PP) * KT / 3! / TT / LF 
1TMP = IT 

IF 1TMP < 80 TIILN TTMP = 80 
I KW = 100 * PP * (FRPM / 1000) A 2.26 / SQR(1TMP) 

FB = FKW / (121 + TT) / LF 
SPR = SP / 14.696 
SPR1 = SRU / KRE 

IF SPR1 > MAXSPR THEN SPR1 = MAXSPR 

IF SPR1 < M1NSPR THEN SPR1 = MINSPR 

IF AUTORE = 1 THEN SPR = SPR1 
SPRU = SPR * 14.696 
COLOR 14, 5 

IF AUTORE = 1 THEN LOCATE 5, 43: PRIN T USING SPRU 


'Temperature Control I »op. 

IF IE = 1 THEN 2950 
IF (TMWL - IT) > 1.5 * MXT THEN 2950 
IF FRPM < 580 11 IEN 2950 
ST = ST1 

IF ST > TT THEN 2530 ELSE 2525 
2525 IF ST < (STP - .04) THEN ST = (STP - .04) 

2530 IF AUTOT = 1 THEN 2540 ELSE 2850 

2540 IF AUTOP = 0 THEN 2610 

IF ABS(PP - SPR) < .15 THEN 2560 ELSE 2570 
2560 IF ABS(TT - S3) < 2! 'll IEN 2610 ELSE 2570 

2570 IF (PP / 13) > .95 * (SPR / S I) THEN 2610 ELSE 2580 

2580 IF IT > ST THEN 2610 ELSE 2590 

2590 ST = TT - .02 

2610 IF (TMWL - ST - MX!) > 0 THEN 2620 ELSE 2650 

2620 ST = TMWL - MXT 
FB = 0 

2650 - IF ST < SAT! II IEN ST = SAT1 

ET = IT - ST 

IF ABS(1T - ST1) < .3 13 IEN ELI = 1 LISE FL1 = 0 
IF ABS(TMWI. - 11) < 24 THEN FI2 = 1 ELSE FL2 = 0 
IF ET < -5 THEN FI 7 = 0 
IF ET > 0 THEN El 7 = I 
FBI 7 = FB * ET * .8 

RIT = RITM1 + KIT * K EGS * DEL * ET 
IF RIT < -FBI 7 THEN RIT = -FBI 7 
IF RIT > (1 - FBI) 11IEN RIT = (1 - FBI 7 ) 

IF (TMWL - 11) > MXT / 2 11IEN 2750 ELSE 2780 
2750 LMT = 1 - (TMWL - IT - MXT / 2) * 2 / MXT 
IF LMT < 1! / LF THEN LMT = 1! / LF 
IF RIT > LMT 11 IEN RIT = LMT 

2780 ALQ = KTGS * (KIT * ET + KDT * (ET - 2 * ETM1 + ETM2) / 2 / DEL) + RIT + FBF 
IF ALQ < 0 11 IEN ALQ = 0 
II 7 ALQ > 1 THEN ALQ = 1 

IF (TMWL - 11) > MXT / 2 11 IEN 2820 ELSE 2830 
2820 IF ALQ > LMT THEN ALQ = LME 
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2830 LCMDS = ALQ * 100 
GOTO 2970 

2850 IF LCMDS > 100 THEN LCMDS = 100 

ALQ = LCMDS / 100 

IF (TMWL - TI) > MXT / 2 THEN 2880 EUSE 2910 
2880 LMT = 1 - (TMWL - IT - MXT / 2) * 2 / MX T 
II ; LMT < I! / LF THEN LMT = 1! / LF 
IE* ALQ > LMT THEN ALQ = LMT 
2910 RIT = ALQ 
GOTO 2970 
2950 ALQ = 0 

LCMDS = 0 

2970 IF ALQ > ALQP + .01 THEN ALQ = ALQP + .01 

DAC(l) = ALQ 
DAC(2) = ALQ 


’Pressure - Reynolds Number Control Ijoop. 

KPGS = 750 * KP / PP / SQR(TI) 

GF1 = 2.725 ♦ CGV * PP / SQR(TI) 

IF PPUSCS > MAXP1 THEN PIE = 1: GOSUB 6000 
IF IE = 1 THEN 3220 

IF ABS(RE - SRE) < .05 THEN FL7 = 1 ELSE FL7 = 0 
IF AUTOP = 1 THEN 3040 EHSE 3180 
3040 EP = PP - SPR 

IF ABS(EP) < .005 THEN El 3 = 1 EHSE FI 3 = 0 
RIP = RIPM1 + EP * KPGS ♦ KIP * DEL 
IF RIP < 0 THEN RIP = 0 
IF RIP > 1! THEN RIP = 1! 

AG VI = KPGS * (KPP * EP + KDP * (EP - 2 * EPM1 + EPM2) / 2 / DEL) + RIP 
IF AGV1 < 0 THEN AG VI = 0 
IF AG VI > 1! THEN AGV1 = 1! 

IF AG VI > .9 THEN AGV2 = AGV2 + .01 
IF AG VI < .7 THEN AGV2 = AGV2 - .01 
IF AGV2 < 0 THEN AGV2 = 0 
IF AGV2 > 1! THEN AGV2 = 1! 

GCMDS = AG VI * 100 
GOTO 3240 

3180 IF GCMDS > 100 THEN GCMDS = 100 

RIP = GCMDS / 100 
AG VI = GCMDS / 100 

IF FRPM < 300 THEN LCMDS = 0 
GOTO 3240 
3220 AG VI = 1! 

GCMDS = UK) 

3240 IF AG VI < AGV1P - .05 THEN AGV1 = AG VIP - .05 

IF AG VI > AG VIP + .05 THEN AGV1 = AGV1P + .05 
DAC(4) = AGV1 * .8 + .2 
DAC(5) = AGV2 * .8 + .2 

’Pan Speed - Mach Number Control I^oop. 

IF IE = 1 THEN 3405 EI.SE 3410 
3405 NCMDS1 = 0 
GOTO 3550 

3410 NCMDS1 = NCMDS 

IF FRPM > MAXN1 THEN NIE = 1: GOSUB 6000 
IF AUTOM = 1 THEN 3422 EHSE 3550 
3422 IF FRPM < 580 THEN AUTOM = 0: GOTO 3435 
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3432 


3435 


3440 


3550 


3560 

3570 

3575 


FKM = 1 - .5 * M 

II- l-'KM < .71 THEN I -KM = .71 
KMC = 520 * SQR(TI) * / IM ’ * 035 

KMGS = KMC * KMM 

IF ALQ > .99 THEN 3432 ELSE 3440 

AUTOM = 0 
NCMDS = FRPM - 500 

IF NCMDS < 0 IHEN NCMDS = 0: BEEP 

COLOR 14, 6 

LOCATE 3, 64: PRINT " MANUAL PRINT 
LOCATE 6, 62: PRINT USING "####."; NCMDS 
GOTO 3550 
EM = SM - M 

IF ABS(EM) < .002 THEN FI A = 1 El.SE FI4 = 0 
RIM = RIMMI + KIM * EM * DEL * KMGS 
IF RIM < -100 THEN RIM = -100 
IF RIM > 100 THEN RIM = 100 
NCMDS! = KPM * EM * KMGS + RIM + FRPM 

IF NCMDS 1 > 5450 THEN NCMDS1 = 5450 
IF NCMDS 1 > MAXN THEN NCMDS1 = MAXN 
EN = NCMDS 1 - FRPM 

IF ABS(EN)<10 THEN FLN=1 ELSE ILN=0 
IF ABS(EN) > 100 THEN 3560 ELSE 3570 
KEN = 5 
GOTO 3575 
KEN = SQR(TT) / 20 
RIN = RINM1 + EN * DEL * KIN 

IF RIN > R1NM1 + KEN THEN RIN = RINM1 + KEN 
IF RIN < RINM1 - KEN THEN RIN = RINMI - KEN 
IF RIN < 0 THEN RIN = 0 
IF RIN > 6000 THEN RIN = 6000 
SRPM = EN * KPN + RIN 

IF SRPM > SRPM1 + 5 THEN SRPM = SRPM1 + 5 
IF SRPM < SRPM1 - 5 THEN SRPM = SRPM1 - 5 
SNRPM = SRPM / 7500 

IF SNRPM > 1! THEN SNRPM = 1! 

IF SNRPM < 0! IHEN SNRPM = 0! 

DAC(6) = SNRPM 


’Fan Sneed Band Warning. 

IF FRPM < 3651 IHEN 3790 ELSE 3830 
3790 IF FRPM > 3549 'IHEN 3800 ELSE 3830 

3800 COLOR 20, 8 

LOCATE 13, 61: PRINT ” SPEED BAND PRINT 
BEEP 


GOTO 3900 
3830 COLOR 14, 6 

LOCATE 13, 61: PRINT H 


PRINT 


’Temperature Gradient Calculation. 

3900 IW = IW + 1 

IF IW = 40 IHEN 3920 ELSE 3980 
3920 IF WLG > 8 THEN WLG = 0 

IF WLG < -8 THEN WLG = 0 
COLOR 14, 3 

LOCATE 20, 30: PRINT USING "###.#"; 12 * WLG 
WLG = 0 
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IW = 0 

.WHO WLG = WI.G + TMWI. - TMWI.l 

’Liquid Hack Pressure Control Ixxip. 

SLQ = SLQSC / 1 4.696 

IF FRPM < <500 THEN SLQ = 9! 

IF FRPM < 400 THEN SIjO = .1 
IF SFQ > 10.2 THEN SLO = 10.2 
SLQSC = SLO * 14.696 
EI.P = PLQ - SLQ 

IF ABS(ELP) < .4 THEN FI -5 = 1 ELSE 11.5 = 0 
IF ELP > .15 WEN ELP = .15 
IF ELP < -.15 THEN ELP - -.15 
RIL = RILM1 + ELP * KIL * DEL 
IF RIL < 0 THEN RIL = 0 
IF RIL > I! WEN RIL = 1! 

ALN = KPI. * ELP + KDL * (ELP - ELPM1) / DEL + RIL 
IF ALN < 0 THEN ALN = 0 
IF ALN > 1! THEN ALN = 1! 

DAC(3) = ALN 
DAC(7) = 1 ! 


’Screen Pressure Drop Warning. 

DLPC = .47 * PP * M * M / (MF) A 6 + .02 

IF DLP > DLPC WEN IT/) = 1 ELSE IT/> = 0 
IF M < .25 WEN FI A = 0 
IF FI/S = 1 WEN BEEP 


’Controller Screen Update. 

JD = JD + 1 

IF JD = 3 WEN 4160 ELSE 4700 
4160 COLOR 14, 3 

LOCATE 21, 30: PRINT USING "###.#"; SAT1 

IF (FL1 * FL2 * AUTOT) = 1 WEN 4190 ELSE 4200 
4190 COLOR 14, 0 

4200 LOCATE 8, 26: PRINT' STRING$(6, 0) 

COLOR 14, 5 

LOCATE 21, 47: PRINT' USING "###.##”; PSUSCS 
IP' IT/) = 1 THEN COLOR 30, 5 
LOCATE 22, 48: PRINT' USING "##.###"; DLP 
IF AUTORE = 0 THEN 4240 
IF (FI .7 * AUTORE) = 1 WEN 4250 ELSE 4260 
4240 IF (FL3 * AUT'OP) = 1 WEN 4150 ELSE 4260 

4250 COLOR 14, 0 

4260 LOCATE 8, 43: PRINT STRING$(6, 0) 

COLOR 14, 2 

LOCATE 5, 13: PRINT!' USING "###.#"; SLQSC 
LOCAT’D 9, 13: PRINT USING "##0"; PLQUSCS 
LOCATE 12, 13: PRINT USING "MO”; (1 - ALN) * 100 
IF FI 25 = 1 THEN 4320 ELSE 4330 
4320 COLOR 14, 0 

4330 LOCATE 8, 13: PRINT’ ST1UNG$(6, 0) 

COLOR 14, 6 

IF (IT A * AUTOM) = 1 THEN 4360 ELSE 4370 
4360 COLOR 14, 0 

4370 LOCATE 8, 62: PRINT' STRING$(6, 0) 
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11 (FI ,N-AUTOM)><) THEN 4374 EI.SE 4376 
4374 COLOR 14,0 

4376 LOCATE 11,62:PRINT STR1N$(6,000) 

COLOR 14, 3 

IF AUTOT = 1 THEN LOCATL 6, 26: PRINT USING "###.# ; SI 
LOCATE 9, 26: PRINT USING "###.#"; TT 
LOCATE 10, 26: PRINT' USING "###.#"; TMWL 
LQ = AI.Q * 100 

LOCATE 12, 26: PRIN T USING "###.#"; LQ 
COLOR 14, 5 

LOCATE 9, 42: PRINT USING "###.##"; PPUSCS 
GV1 = AG VI * 100 

LOCATE 12, 43: PRIN T USING "###.#"; GV1 
GV2 = AGV2 * 100 

LOCATE 13, 43: PRINT USING "###.#”; GV2 
LOCATE 10, 42: PRINT USING "###.##"; RE 
COLOR 14, 6 

LOCATE 9, 62: PRINT USING "#.###"; M 
LOCATE 10, 62: PRINT USING ”####."; FRPM 
LOCATE 12, 62: PRIN T USING SNRPM * 100 


JD = 0 


4700 GOSUH 70(X) 

’Program Timing Counter. 

K = K + 1 

IF K = 600 THEN HELP LINE 4800 

K = 0 


’Selling Previous Cycle Values. 

4800 AG VIP = AG VI 
ALQP = ALQ 
ELPM1 = LLP 
EPM2 = EPM1 
EPMI = EP 
Cra2 = L'TMl 
ETM1 = ET 
RILM1 = RIL 
RIMM1 = RIM 
RINM1 = RIN 
RIPM1 = RIP 
RITM1 = RIT 
SRPM1 = SRPM 
SIT = ST 
GOTO 1000 


’Conlrollcr Input Subroutine. 

S(HX) AS = INKLYS 

IF AS = "" IHEN 5‘W9 


5008 

5009 

5010 


IF' AS = "Q" OR AS = ”q" THEN 7500 

IF AS = "E" OR AS = "e" THEN 6000 

IF AS = -D" OR AS = ”d" THEN 5950 

IF' XX > 5 H1EN XX = 5 

IF AS = CI1R$(13) THEN 5008 EINE 5009 

IF XX = 5 THEN 5600 LINE 5999 

IF J > 0 THEN 5010 LINE 5014 

IF AS = "0" OR AS = "I” OR AS = "2" OR AS = M 3" 


OR AS = "4" OR AS = ”5" 


OR AS 
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= "6" OR AS = "7* OR AS = "8" OR AS = "9" OR AS = THEN 5100 ELSE 5999 
5014 COLOR 16, 3 

IF AS = T M OR AS = "l” THUN J = 1: LOCATE 15, 36: PRINT M ,K M : PRINT ; 

IF AS = "L" OR AS = "l" THEN J = 5: LOCATE 16, 36: PRINT \opn": PRINT ; 
COLOR 16, 5 

IF AS = "P M OR AS = "p M THEN J = 2: LOCA1T 15, 54: PRINT ",psia": PRINT; 

IF AS = "R" OR AS = V THEN J = 4: LOCATE 16, 54: PRINT ",miln": PRINT; 

IF AS = "G" OR A$ = "g" THEN J = 6: LOCATE 17, 54: PRINT *\opn": PRINT; 

IF AS = "C" OR AS = M c" THEN J = 8: LOCATE 18, 54: PRINT H ,in H : PRINT ; 
COLOR 16, 6 

IF AS = "M" OR AS = "m" THEN J = 3: LOCATE 15, 71: PRINT M ,Mach M : PRINT; 
IF AS = "N" OR AS = "n" THEN J = 7: LOCA1T 16, 71: PRINT \rprn": PRINT; 
COLOR 16, 2 

IF AS = M ir OR AS = ”b" THEN J = 9: LOCATE 15, 19: PRINT ",psi M : PRINT ; 

IF J = 1 THEN LU = 15: MU = 30 

IF J = 2 ITIEN LU = 15: MU = 48 

IF J = 3 THEN LU = 15: MU = 65 

IF J = 4 11IEN LU = 16: MU = 48 

IF J = 5 THEN LU = 16: MU = 30 

IF J = 6 THEN LU = 17: MU = 48 

IF J = 7 THEN LU = 16: MU = 65 

IF J = 8 THEN LU = 18: MU = 48 

IF J * 9 THEN LU = 15: MU = 13 

IF J > 0 11IEN 1IEEP 

GOTO 5999 

5100 IF J = 1 THEN 5105 

IF J = 2 THEN 5405 
IF J = 3 THUN 5205 
IF J = 4 THEN 5405 
IF J = 5 THEN 5405 
IF J = 6 THEN 5405 
IF J = 7 THEN 5350 
IF J = 8 THEN 5405 
IF J = 9 THEN 5105 

’Assemble Temperature and Liquid Hack Pressure Inputs. 

5105 XX = XX + 1 

D1 = 1: D2 = 2: D3 = 3: D4 = 4: D5 = 5 
GOSUH 5500 

AA = III * 100 + 1 12 + 10 + 1 13 + 1 15 / 10 
GOTO 5999 

’Assemble Mach Number Input. 

5205 XX = XX + 1 

D1 = 1: D2 = 3: D3 = 4: D4 = 2: D5 = 5 
GOSUH 5500 

AA = 112 / 10 + 113 / 100 + 115 / 1000 
GOTO 5999 

’Assemble Fan Speed Input. 

5350 XX = XX + 1 

D1 = 1: D2 = 2: D3 = 3: D4 - 5: D5 = 4 
GOSUH 5500 

AA = III * KMX) + 112 * 100 + 1 13 * 10 + 115 
GOTO 5999 
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’Assemble Pressure, Reynolds Number, Injection Valve, Exhaust Valve, and Chord Inputs. 

5405 XX = XX + 1 

D1 = 1: D2 = 2: D3 = 4: 1)4 = 3: D5 = 5 
GOSUB 5500 

AA = III * 10 + 112 + 113 / 10 + 115 / 100 
GOTO 5999 


55(H) 


IE 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

RETURN 


J = 1 THEN COLOR 0, 3 
J = 2 THEN COLOR 0, 5 
J = 3 THEN COLOR 0, 6 
J = 4 THEN COLOR 0, 5 
J = 5 THEN COLOR 0, 3 
J = 6 THEN COLOR 0, 5 
J = 7 THEN COLOR 0, 6 
J = 8 THEN COLOR 0, 5 
J = 9 THEN COLOR 0, 2 
XX = 1)1 THEN 111 = VAIXAS) 
XX = 1)2 THEN 1 12 = VAI<A$) 
XX = 1)3 THEN 113 = VAUA$) 


LOCATE LU, MU + 

LOCATE LU, MU + 

LOCATE LU, MU + 

XX = 1)4 THEN LOCATE LU, MU + 1)4: PRINT CIIR$(46); 

XX = D5 THEN 115 = VAI<AS>: LOCATE LU, MU + D5: PRINT USING "#";II5 


Dl: PRINT USING "# H ;TI1 
D2: PRINT USING "#”;II2 
1)3: PRINT USING "#";II3 


’Set Temperature Input. 

5600 IF J = 1 THEN 5620 ELSE 5660 

5620 ST1 = AA 

IF ST1 > MAXT THEN ST1 = MAXT 
IF ST1 < MINT THEN ST1 = MINT 
IF ST1 < SAT1 'IT1EN ST1 = SAT1 
SIT = IT 
COLOR 14, 3 

LOCATE 3, 29: PRINT " AUTO PRINT ; 
A LO OT = 1 

LOCATE 5, 26: PRINT USING "###.#"; ST1 


’Set Pressure Input. 

5660 IF J = 2 'THEN 5670 ELSE 5730 

5670 SP = AA 

IF SP < M1NP THEN SP = MINP 
IF SP > MAXP THEN SP = MAXP 
COLOR 14, 5 

LOCATE 3, 46: PRIN T " AUTOP PRINT ; 

LOCATE 6, 43: PRIN'I : PRIN T ; 

LOCATE 5, 43: PRIN T USING "##.##"; SP 
AUTOP = I 
AUTORE = 0 


’Set Mach Number Input. 

5730 IF J = 3 THEN 5740 ELSE 5750 

5740 SM = AA 

IF SM < M1NM THEN SM = MINM 
IF SM > MAXM THEN SM = MAXM 
COLOR 14, 6 

I OCATE 3, 64: PRIN T " AUTO ": PRIN'I’ ; 

LOCATE 6, 62: PRIN'I 

AUTOM = 1 
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LOCATE 5, 62: PRINT USING "#.###"; SM 


’Set Reynolds Number Input. 

5750 IF J = 4 THEN 5760 ELSE 5770 

5760 SRE = AA 

IF SRE > MAXRE T HEN SRE = MAXRE 
IF SRE < MINRE THEN SRE = MINRE 
COLOR 14, 5 
AUTORE = 1 
AUTOP = 1 

LOCATE 6, 43: PRINT USING "##.##"; SRE 
LOCATE 3, 46: PRINT " AUTORE ": PRINT ; 


’Set Liquid Injection Valve Input. 

5770 IF J = 5 THEN 5780 ELSE 5810 

5780 LCMDS = AA 
COLOR 14, 3 

LOCATE 3, 29: PRINT " MANUAL ": PRINT ; 

LOCATE 5, 26: PRIN1 

LOCATE 6, 26: PRIN’I 

AUTOT = 0 

’Set Gas Exhaust Valve Input. 

5810 IF J = 6 THEN 5820 ELSE 58(4) 

5820 GCMDS = AA 
COLOR 14, 5 

LOCATE 3, 46: PRINT " MANUAL ": PRINT ; 

AUTOP = 0 
AUTORE = 0 

LOCATE 5, 43: PRIN1 : PRIN T ; 

LOCATE 6, 43: PRIN1 : PRINT ; 

’Set I'an Speed Input. 

5860 IF J = 7 THEN 5870 ELSE 5935 

5870 NCMDS = AA 

IF NCMDS > MAXN THEN NCMDS = MAXN 
COLOR 14, 6 

LOCATE 5, 62: PRIN1 

LOCATE 3, 64: PRINT " MANUAL PRIN T ; 

LOCATE 6, 62: PRINT USING "####."; NCMDS 
AUTOM = 0 

’Set Chord Ixngth Input. 

5935 IF J = 8 THEN 5937 ELSE 5943 

5937 CHIN = AA 

IF CHIN < MINCHIN THEN CHIN = MINCIIIN 
IF CHIN > MAXCIIIN THEN CHIN = MAXCIIIN 
COLOR 14, 5 

LOCATE 20, 48: PRINT USING "##.##"; CHIN 
Cl I = CHIN * .0254 

’Set Liquid Back Pressure Input. 

5943 IF J = 9 THEN 5944 ELSE 5950 

5944 SLQSC = AA 

IF SLQSC > MAXLQSC THEN SLQSC = MAXLQSC 
COLOR 14, 2 

LOCATE 5, 13: PRINT USING "###.#"; SLQSC 
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5950 IF J = 0 THEN 5962 

IF J = 1 THEN COLOR 14, 3 
IF J = 2 TIIliN COLOR 14, 5 
IF J = 3 TIIliN COLOR 14, 6 
IF J = 4 THEN COLOR 14, 5 
IF J = 5 TIIliN COLOR 14, 3 
IF J = 6 TIIliN COLOR 14, 5 
IF J = 7 TIIliN COLOR 14, 6 
IF J = 8 THEN COLOR 14, 5 
IF J = 9 THEN COLOR 14, 2 
LOCATE LU, MU + 1: PRINT " PRINT ; 

5962 IF J > 0 THEN BEEP 

IF IE = 1 THEN COLOR 0, 6 
IF IE = 1 THEN BEEP: LOCATE 20, 60: PRINT " 

IF PIE = 1 THEN LOCATE 21, 60: PRIN T " 

IF NIE = 1 THEN LOCATE 21, 60: PRINT - 

J = 0 
XX = 0 
IE = 0 

5999 RETURN 

’Emergency Stop Subroutine. 

6000 IE = 1 
COLOR 20, 8 

LOCATE 20, 60: PRINT " EMERGENCY STOP ": PRINT ; 

COLOR 4, 8 

IF PIE = 1 11IEN LOCATE 21, 60: PRIN T " Pressure Limit 

IF NIE = 1 THEN LOCATE 21, 60: PRINT " Fan RPM Limit " 

COLOR 4, 8 
AUTOT = 0 
AUTOP = 0 
AUTOM = 0 
AUTORE = 0 

LOCATE 3, 29: PRINT " MANUAL M : PRINT ; 

LOCATE 3, 46: PRINT " MANUAL *: PRINT ; 

LOCATE 3, 64: PRINT " MANUAL ": PRINT ; 

RETURN 


PRINT ; 

PIE = 0 
B : NIE = 0 


’Simulation of 0.3-m TCT Dynamics. 

7000 WG = 4375 * PP / TT 
WT = 3200! 

CM = (5.5 * TT - .008 * TV A 2) / 1000 
CP = 1.04 


CV = .75 

TMC = .28 * TT A 1.2 / PP A .7 / M * .7 


LMFL = ALQ * LF 

GMFT. = (AG VI + AGV2) * CGV * 2.725 * PP / SQR(TT) _ 

IF PP < 1.5 THEN GMFL = (AGV1 + AGV2) * PP / SQR(TT) 21.8 (2 - (1.5 / PP) 1.7) 

WLHT = (TMWL - rn * WT * CM / TMC 

HEAT = WLIIT + FKW - LMFL * (121 + CV * TO - GMFL * (CP - CV) * TT 
DPP = DEL * PP * HEAT / TT / WG / CV + DEL * (LMFL - GMFI.) * PP / WG 


PPN = PP + DPP 

IF PPN > 6! 'THEN PPN = 6! 
IF PPN < 1! THEN PPN = 1! 


47 



TMWLN = (D1ZL / (2 * TMC + DEL)) * (IT + TIM) + ((2 * TMC - DEL) / (2 * TMC + DEL)) * 
TMWL 

DFR1 = (1 / (DEL A 2 + LI 2 * DEL + .8)) 

DFR2 = (DEI. A 2) * (SRPM + SRPM1 * 2 + SKPM2) - FRPM * (2 * DEL A 2 - 1.6) 

DFR3 = DFR2 - FRPM1 * (DEI. * 2 - 1.12 * DEL + .8) 

FRPMN = DERI * DFR3 

DMT = DEL * HEAT / WG / CV 

TYN = TT + DTT 

ALF1 = (ALN - .2) * 5 

PIX) = 7.3 + ALF1 * 3 + ALO * .7 

PLQN = DEL / (2 + DEL) * (PLO + PLOl) + (2 - DEL) / (2 + DEL) * PLQ 
IF PLQN < 7.3 MMIEN PLQN = 7.3 
PP1 = PP 
r IT1 = 1M 1 
TMWL1 = TMWL 
VRVMX = ETPM 
SRPM2 = SRPM1 
PLOl = PLO 
IT = rm 

IF IT < SAT I THEN TP = SAT1 
PP = PPN 
TMWL = TMWLN 
FRPM = FRPMN 
PLQ = PLQN 
RETURN 

7500 CLS 

SYSTEM 
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Appendix C 
Source Code Variables 


Limits Input Keyboard Commands 

A, a Accepts the tunnel maximum-minimum limits and starts the 0.3-m TCT controller. 

D,d Delete the previous input keys shown on the screen and not yet executed. 

M, m Input maximum Mach number limit. Maximum Mach number limit format is #.###, 

range to 0.15 to 0.995. 

N, n Input maximum fan speed limit. Maximum fan speed limit format is ####., range 

is 0 to 5600 rpm. 

P,p Input maximum tunnel total pressure limit. Maximum tunnel total pressure limit 
format is ##.##, range is 14.7 to 88.0 psia. 

R, r Input maximum Reynolds number limit. Maximum Reynolds number limit format is 

##.##, range is 1 to 50 million. 

S, s Accepts the user set maximum-minimum limits and starts the 0.3-m TCT controller. 

T, t Input minimum temperature limit. Minimum temperature limit format is ###.#, 

range is 77 to 340 K. 


Controller Input Keyboard Commands 

B, b Input LN 2 back pressure set point. LN 2 back pressure set point format is ###.#, 

range is LN 2 pressure to 150 psia. 

C, c Input mean aerodynamic chord. Chord format is ##.##, range is 0.39 to 15.75 

inches. 

D, d Delete the previous input keys memorized on the screen and not yet executed. 

G,g Input GN 2 discharge valve area. Takes the pressure controller and Reynolds number 
controller to manual control mode. GN 2 discharge valve area format is ##.##, range 
is 99.99%= full open to 0%=closed. 

LJ Input LN 2 injection valve area. Takes the temperature controller to manual control 
mode. LN 2 injection valve area format is ##.##, range is 99.99%=full open to 
0%=closed. 
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M, m Input Mach number set point. Takes fan speed to automatic Mach number control 

mode. Mach number set point format is #.###, range is 0.150 to 0.995. 

N, n Input fan speed set point. Takes the fan speed to manual Mach number control 

mode. Fan speed set point format is ####., range is 0 to 5600 rpm. 

P,p Input tunnel total pressure set point. Takes the controller to automatic pressure 

control mode and manual Reynolds number control mode. Total pressure set point 
format is ##.##, range is 14.7 to 88 psia. 

R,r Input Reynolds number set point. Takes the pressure controller to automatic 
Reynolds number control mode by generating the required pressure set point. 
Reynolds number set point format is ##.##, range is 1 to 50 million. 

T,t Input temperature set point. Takes the controller to automatic temperature control 

mode. Temperature set point format is ###.#, range is saturation temperature to 340 
K. 


Program Variables 


AA Keyboard input function. Combines the integer inputs HI, H2, H3, and H5 

with proper decimal scaling for use in program calculation. 

AEMP Temporary data register output in digital to analog conversion routine to 

check status. 


AG VI 

AGV1P 

AGV2 

ALF1 

ALN 

ALO 


Area of GN, discharge valve #1. Full open AGV1=1 and closed AGV1=0. 
AG V 1 =KPGS*(KPP*EP+KDP*(EP-2*EPMl+EPM2/2/DEL)+RIP 

Area of GN, discharge valve #1, previous cycle value. 

Area of GN 2 discharge valve #2. Full open AGV2=1 and closed AGV2=0. 
Valve #2 opens only when AGV1 > 90% open and starts closing when 
AGV1 < 70% open. Valve #2 moves at a rate of 1% per cycle. 

Simulator variable. Variable related to the area of LN 2 back pressure valve 
used in calculating the estimated LN 2 back pressure. 

ALFl=(ALN-.2)*5 

Area of LN 2 back pressure control valve. Full open ALN=1 and closed 
ALN=0. ALN=KPL*ELP+KDL*(ELP-ELPM1)*DEL+RIL 

Area of LN 2 injection valve. Full open ALQ=1 and closed ALQ=0. 

AIX3=KTGS*(KPT*F.T+KDT*(RT-2*ETMl+F.TM2)/2/DEL)+RIT+FBF 


ALQP Area of LN 2 injection valve, previous cycle value. 
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AUTOM 

AUTOP 

AUTORE 

AUTOT 

A$ 

BEMP 

BSD 

B$ 

CEMP 

C$ 

CGV 

CH 

CHIN 

CLOV 

CM 

CP 

CRPM 


Mach number controller: 

AUTOM=l automatic Mach number control. 

AUTOM=0 manual Mach number control. 

Pressure controller: 

AUTOP=l automatic pressure control. 

AUTOP=0 manual pressure control. 

Reynolds number controller: 

AUTORE= 1 automatic Reynolds number control. 

AUTORE=0 manual Reynolds number control. 

Temperature controller: 

AUTOT=l automatic temperature control. 

AUTOT=0 manual temperature control. 

Input variable from the keyboard buffer. 

Temporary status register output in digital to analog conversion routine to 
check status. 

Dummy calculation to create a time delay between resetting the digital to 
analog conversion and reading its status. Used only during initialization of 
digital to analog conversion. 

Input variable from the keyboard buffer. 

Temporary status register output in digital to analog conversion subroutine. 
Input variable from the keyboard buffer. 

Coefficient for GN 2 exhaust valve. 

Aerodynamic chord in meters. 

Aerodynamic chord in inches. Loaded through the keyboard using the "C" 
command. Default value is CHIN =0.1800*0.0254. 

Flow coefficient of LN 2 injection valve. 

Simulator variable. Specific heat of tunnel metal. 
CM=(5.5*TT-0.008*TT A 2)/1000 

Simulator variable. Specific heat of nitrogen at constant pressure. 
Simulator variable. Fan speed in rpm with a minimum limit of 50 rpm. 
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CRM 

CV 

DAC(I) 

DAC(2) 

DAC(3) 

DAC(4) 

DAC(5) 

DAC(6) 

DAC<7) 

D1 

D2 

D3 

D4 

D5 

DEL 

DFR1 

DFR2 

DFR3 

DLP 


Simulator variable. Variable used in calculating Mach number as a function 
of fan speed, pressure, and temperature. 
CRM=2.77-CRPM/157/SQR(TT)*PP A 0.035 

Simulator variable. Specific heat of nitrogen at constant volume. 

Control variable ALQ output. 

Control variable ALQ output. 

Control variable ALN output. 

Control variable AGV1 output. 

Control variable AGV2 output. 

Control variable SNRPM output. 

RPM/MACH controller output for rheostat control. DAC(7)=1 for normal 
operation. DAC(7)=0 for rheostat control. 

Integer used to print properly the variable being loaded through the keyboard. 

Integer used to print properly the variable being loaded through the keyboard. 

Integer used to print properly the variable being loaded through the keyboard. 

Integer used to print properly the variable being loaded through the keyboard. 

Integer used to print properly the variable being loaded through the keyboard. 

Cycle time step. 

Simulator variable. Variable used in calculating the estimated fan speed rpm. 
DFR 1 =( 1 /(DEL A 2+ 1 . 1 2*DEL+0.8)) 

Simulator variable. Variable used in calculating the estimated fan speed rpm. 

DFR2=(DEL A 2)*(SRPM+SRPM1*2+SRPM2)-FRPM*(2*DEL A 2-1.6) 

Simulator variable. Variable used in calculating the estimated fan speed rpm. 
DFR3=DFR2-FRPM 1 *(DEL A 2-1 . 1 2*DEL+.8) 

Difference in pressure across the screens of the tunnel settling chamber, 

(psia). 
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DLPC 

DPP 

DRPM 

DR1 

DR2 

DTT 

ELP 

ELPM1 

EM 

EN 

EP 

EPM1 

EPM2 

ET 

ETM1 

ETM2 

E(l) 

E(2) 

E(3) 

E(4) 

E(5) 


Maximum sale screen pressure drop allowed. 
DLPC=0.47 *PP*M*M/(MF) a 6+0.02 

Simulator variable. Estimated change in tunnel total pressure for one cycle. 

DPP=DEL*PP*IIEAT/TT/WG/CV+DEL*(LMFL-GMFL)*PP/WG 

Fan speed in rpm with a minimum limit of 100 rpm. 

Command output high integer buffer in digital to analog conversion routine, 
(12 bit). 

Command output low integer buffer in digital to analog conversion routine, 
(4-bit). 

Simulator variable. Estimated change in tunnel temperature for one cycle. 
DTT=DEL*HEAT/WG/CV 

Error in LN 2 back pressure control, (atm). ELP=SLQ-PLQ 
Error in LN 2 back pressure control, previous cycle value. 

Error in Mach number control. EM=SM-M 
Error in fan speed control, (rpm). EN=NCMDS1-FRPM 
Error in total pressure control, (atm). EP=PP-SPR 
Error in total pressure control, previous cycle value. 

Error in total pressure control, twice previous cycle value. 

Error in temperature control, (kelvin). ET=TT-ST 
Error in temperature control, previous cycle value. 

Error in temperature control, twice previous cycle value. 

Input millivolts from total pressure sensor. 

Input millivolts from static pressure sensor. 

Input millivolts from total temperature thermocouple. 

Input millivolts from the tunnel metal wall temperature thermocouple. 

Input millivolts from fan speed sensor. 
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m 

E(7) 

Eli 

E12 

FB 

FBF 

FF 

FKM 

FKW 

FLO 

FL1 

FL2 

FL3 

FL4 

FL5 


Input millivolts from LN 2 pressure transducer. 

Input millivolts from the screen pressure drop transducers. 

Output from analog to digital conversion of the tunnel total pressure for 
monitor display. 

Output from analog to digital conversion of the tunnel static pressure for 
monitor display. 

Fan bias, an equivalent of FKW in LN 2 flow. FB=FKW/(121+TT)/LF 
Product of FB and FF. FBF=FB*FF*0.8 

Feed forward logic integer. FF=1 the fan bias is fed forward as LN 2 . FF=0 
the fan bias is not fed forward as LN 2 . The feed forward is off during 
tunnel warm-ups. 

Variable used in calculating the tunnel fan speed/test section coupling. 
FKM=l-0.5*M 


Estimated fan power released to gas, (kilowatts). 
FKW=100*PP*(FRPM/1000 ) a 2.26/SQR(TTMP) 

Logic flag. FLO-1 when the set limit maximum equals the tunnel minimum 
value. This condition is an unacceptable condition to run the tunnel 
controller. FL0=0 corresponds to acceptable tunnel set limits to run the 
tunnel controller. 

Logic flag. FL1 = 1 when the error in the temperature control is less than 0.3 
K. FL1=0 corresponds to a larger temperature control error. 

Logic flag. FL2=1 when metal temperature is within 24 K of the tunnel 
temperature set point. FL2=0 corresponds to a larger temperature difference. 
FL1*FI2=1 corresponds to the temperature of the tunnel at the temperature 
set point. 

Logic flag. FL3=1 when the error in pressure control is less than 0.005 atm. 
FL3=0 corresponds to a larger pressure control error. 

Logic flag. FL4=1 when the error in Mach number is less than 0.002. 

FL4=0 corresponds to a larger Mach number control error. 

Logic flag. FL5=1 when the error in the LN 2 back pressure control is less 
than 0.4 atm. FL5=0 corresponds to a larger LN 2 back pressure error. 
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FL6 


Logic flag. FL6=1 when the pressure difference across the tunnel screens is 
unsafe. FL6=0 corresponds to a safe pressure difference. 

FL7 Logic flag. FL7=1 when the error in the Reynolds number control is less 

than 0.05. FL7=0 corresponds to a larger Reynolds number control error. 

FL10 Logic flag. FL10=1 when a sensor fails. FL10=0 corresponds to normal 

operation. A sensor failure will cause an emergency stop. 

FRPM Fan speed in rpm. 

FRPMN Simulator variable. Estimated fan speed. FRPMN=DFR1*DFR3 

GCMDS GN 2 discharge valve area. Loaded through the keyboard using the "G" 
command. 

GFI Mass flow from the tunnel. GF1=2.725*CGV*PP/SQR(TT) 

GMFL Simulator variable. GN 2 exhaust mass flow. 

GV1 Product of AG VI and 100. GV1=AGV1*100 

GV2 Product of AGV2 and 100. GV2=AGV2*100 

HI First integer of an input keyed during a variable input. 

H2 Second integer of an input keyed during a variable input. 

H3 Third integer of an input keyed during a variable input. 

H5 Fourth integer of an input keyed during a variable input. 

HEAT Simulator variable. The total heat flow of the tunnel. 

HEAT=WLHT+FKW-LMFL*(121+CV*TT)-GMFL*(CP-CV)*TT 

HI High byte to digital to analog conversion buffer of the digital to analog 

conversion routine. 

I Integer cycle counter used to display the controller screen borders. 

IE Logic flag. IE=1 corresponds to an emergency stop. IE=0 corresponds to 

normal operation. 

IK Integer cycle counter used for output of digital to analog conversion 

variables. 

IW Integer cycle counter for integrating the wall temperature gradient. 
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J 



JD 

JJ 

K 

KDL 

KDP 

KDT 

KEN 

KIL 

KIM 

KIN 

KIP 

KIT 

KMC 

KMGS 

KMM 

KP 

KPGS 

KPL 

KPM 

KPN 

KPP 

KPT 


Integer variable that identifies the input variables. 

Integer cycle counter for displaying certain variables on the screen. 

Integer cycle counter for calculation of BSD. 

Simulator variable. Integer cycle counter for a 600 cycle beep. 

LN 2 back pressure control derivative gain. 

Pressure control derivative gain. 

Temperature control derivative gain. 

A gradient term limiting the rpm of the fan speed. 

LN 2 back pressure control integral gain. 

Mach number control integral gain. 

Fan speed control integral gain. 

Pressure control integral gain. 

Temperature control integral gain. 

Constant corresponding to tunnel fan speed/test section coupling in the Mach 
number controller. KMO=520*SQR(1T)*FKM/PP A 0.035 

Mach number control gain schedule function. KMGS=KMC*KMM 

Mach number control gain. 

Pressure control gain. 

Pressure control gain schedule function. KPGS=750*KP/PP/SQR(TT) 

LN 2 back pressure control proportional gain. KPL=0.2 
Mach number control proportional gain. 

Fan speed control proportional gain. 

Pressure control proportional gain. 

Temperature control proportional gain. 
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KRE 

KT 

KTGS 

LC 

LCMDS 

LDPO 

LF 

LMFL 

LMT 

LO 

LO 

LR 

LU 

M 

MAXCH 

MAXCHIN 

MAXLQSC 

MAXM 

MAXM1 

MAXN 


Constant used for evaluating Reynolds number related functions. 

KRE=637 14*CH*M/TT A 1 ,4/*(MF) A 2. 1 

Temperature control gain factor. 

Temperature control gain schedule function. 
KTGS=DRPM*SQR(PP)*KT/3.0/TT/LF 

Screen column number used in LOCATE statements. 

LN 2 valve area. Loaded through the keyboard using the "L" command. 

Difference between LN 2 back pressure and tunnel pressure. 

LDPQ=PLQ-PP 

LN 2 flow when LN 2 injection valve is full open. 
LF=CLQV*0.8676*SQR(LDPQ) 

Simulator variable. LN 2 injected mass flow. LMFL=ALQ*LF 

Limiter for LN 2 flow. LMT= 1 <TMWL-TT-MXT/2)*2/MXT 

An output driving the digital to analog conversion low buffer for final 
command. Low byte data register in digital to analog conversion routine. 

Product of ALQ and 100. LQ=ALQ*100 

Screen row number used in LOCATE statements. 

Screen column number used in LOCATE statements. 

Tunnel flow Mach number in the test section based on the difference 
between total and static pressure. 

Maximum aerodynamic chord length limit in meters. 

Maximum aerodynamic chord length limit in inches. 

Tunnel maximum LN 2 back pressure. 

User set maximum Mach number limit. Must be less than MAXM1. 
Tunnel maximum Mach number limit. 

User set maximum fan speed limit. Must be less than MAXN1. 
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MAXN1 

MAXP 

MAXP1 

MAXRE 

MAXRE1 

MAXSPR1 

MAXT 

MAXT1 

MF 

MINCH 

MINCHIN 

MINM 

MINM1 

MINN 

MINN1 

MINP 

MINP1 

MINRE 

MINRE1 

MINSPR1 

MINT 

MINT1 


Tunnel maximum fan speed limit. 

User set maximum pressure limit. Must be less than MAXP1. 

Tunnel maximum pressure limit. 

User set maximum Reynolds number limit. Must be less than MAXRE1. 
Tunnel maximum Reynolds number limit. 

Maximum pressure set point estimated from the Reynolds number set point. 
Valid in automatic Reynolds number control. MAXSPR1=MAXP1/14.696 

Maximum temperature limit. 

Tunnel maximum temperature limit. 

Isentropic function of the tunnel flow Mach number. MF=(1+0.2*M*M) 
Minimum aerodynamic chord length in meters. 

Minimum aerodynamic chord length in inches. 

Minimum Mach number limit. 

Tunnel minimum Mach number limit. 

Minimum fan speed limit. 

Tunnel minimum fan speed limit. 

Minimum pressure limit. 

Tunnel minimum pressure limit. 

Minimum Reynolds number limit. 

Tunnel minimum Reynolds number limit. 

Minimum pressure set point estimated from the Reynolds number set point. 
Valid in automatic Reynolds number control mode. 

MINSPRUMINP 1/1 4.696 

User set minimum temperature limit. Must be greater than MINT1. 

Tunnel minimum temperature limit. 
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MU Screen column number used in the LOCATE statements. 

MXT Maximum safe temperature difference allowed between tunnel gas and tunnel 

metal wall. 

NCMDS Fan speed set point. Loaded through the keyboard using the "N" command. 

NCMDS1 Fan speed set point derived from NCMDS or other safety fan speed 

commands. 

NIE Logic flag. NIE=1 when the tunnel exceeds the fan speed rpm limit. 

NIE=1 will causes an emergency stop of the tunnel. NIE=0 corresponds to 
normal operation. 

PIE Logic flag. PIE=1 when the tunnel exceeds the total pressure limit. PIE=1 

will cause an emergency stop of the tunnel. PIE=0 corresponds to normal 
operation. 

PLO Simulator variable. Variable used in calculating the estimated LN 2 back 

pressure. PL0=7.3+ALF1 *3+ALQ*.7 

PLOl Simulator variable. Variable used in calculating the estimated LN 2 back 

pressure, previous cycle value. 

PLO LN 2 back pressure, (atm). 

PLQN Simulator variable. Estimated LN 2 back pressure. 

PLQN=DEL/(2+DEL)*(PL0+PL01)+(2-DEL)/(2+DEL)*PLQ 

PLQUSCS LN 2 back pressure, (psia). PLQUSCS=PLQ* 14.696 

PP Tunnel total pressure, (atm). 

PPN Simulator variable. Estimated tunnel total pressure. PPN=PP+DPP 

PP1 Simulator variable. Tunnel total pressure, previous cycle value. 

PPUSCS Tunnel total pressure, (psia). PPUSCS=PP* 14.696 

PS Static pressure, (atm). 

PSUSCS Static pressure, (psia). PSUSCS=PS* 14.696 

RE Flow Reynolds number based on aerodynamic chord. RE=KRE*PP 

READ.WAIT Digital to analog software read command. 
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RIL 

RILM1 

RIM 

RIMM1 

RIN 

RINM1 

RIP 

RIPMI 

RIT 

RITM1 

SAT 

SAT1 

SLO 

SLQSC 

SM 

SNRPM 

SP 

SPR 

SPR1 


LN 2 back pressure control integral error, (atm-sec). 

RIL=RILM 1 +ELP*KIL*DEL 

LN 2 back pressure control integral error, previous cycle value. 

Mach number control integral error, (Mach-sec). 
RIM=RIMM1+KIM*EM*DEL*KMGS 

Mach number control integral error, previous cycle value. 

Fan speed control integral error, (rpm-sec). RIN=RINM1+EN*DEL*KIN 

Fan speed control integral error, previous cycle value. 

Tunnel total pressure control integral error, (atm-sec). 

RIP=RIPM I +EP*KPGS*KIP*DEL 

Tunnel total pressure control integral error, previous cycle value. 

Temperature control integral error, (K-sec). 

RIT^RITM 1 +KIT*KTGS*DEL*ET 

Temperature control integral error, previous cycle value. 

Nitrogen saturation temperature based on Jacobsens data. 
SAT=50+27.34*PS A 0.296 

Nitrogen saturation temperature applied to tunnel static pressure. 
SAT1=SAT*MF 

LN 2 back pressure control set point, (atm). SLQ=SLQSC/14.696 

LN 2 back pressure set point. Loaded through the keyboard using the "B" 
command. 

Mach number set point. Loaded through the keyboard using the "M" 
command. 

Fan speed rheostat drive command normalized to one. SNRPM=SRPM/7500 

Tunnel total pressure set point. Loaded through the keyboard using the "P" 
command. 

Tunnel total pressure set point, (atm). SPR=SP/14.696 

Tunnel total pressure set point estimated from the Reynolds number set point. 
Valid in automatic Reynolds number control mode. SPR1=SRE/KRE 
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SPRU 


Tunnel total pressure set point estimated from the Reynolds number set point, 
(psia). SPRU=SPR* 14.696 

SRE Reynolds number set point. Loaded through the keyboard using the "R" 

command. 

SRPM Fan speed command from control law. SRPM=EN*KPN+RIN 

SRPM1 Fan speed command from control law, previous cycle value. 

ST Tunnel use temperature set point. Based on a safe temperature for given 

tunnel conditions. Also for conditions when both TT and PP are high and 
the gas mass in the tunnel is inadequate to reach the final temperature set 
point. 

ST1 Tunnel temperature set point. Loaded through the keyboard using the "T" 

command, (kelvin). 

STP Tunnel use temperature set point, previous cycle value. 

TMC Simulator variable. Metal time constant of heat release. 

TMC=948/TT A . 12/PP\8/M A .8 

TMWL Tunnel metal wall temperature, (kelvin). 

TMWLN Simulator variable. Estimated tunnel metal wall temperature. 

TMWl>J=(DRiy(2*TMC+DF.I.))*(TT+lTl)+((2*TMC-DP.I.)/(2*TMC+DF.L))*TMWL 

TMWL1 Tunnel metal wall temperature, previous cycle value. 

TT Temperature of tunnel gas, (kelvin). 

TTN Simulator variable. Estimated tunnel gas temperature. TTN=TT+DTT 

TT1 Simulator variable. Temperature of tunnel gas, previous cycle value. 

TTMP Temperature of tunnel gas with a minimum limit of 80 K. 

WG Simulator variable. Mass of nitrogen gas in the tunnel. WG=4375*PP/TT 

WLG Rate of change for the tunnel metal wall temperature. Wall temperature 

gradient. WLG=WLG+TMWL-TMWL1 

WLHT Simulator variable. Rate of heat release from the metal walls in the tunnel. 

WLHT=(TMWL-TT)*WT*CM/TMC 

WRITE. WAIT Digital to analog software write command. 
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WT 


Simulator variable. Mass of the tunnel metal exposed to the tunnel flow. 
XDLP Screen pressure drop transducer sensitivity. 

XFRPM Fan speed sensor sensitivity. 

XPLQ LN 2 back pressure sensor sensitivity. 

XPP Tunnel total pressure transducer sensitivity. 

XPS Static pressure transducer sensitivity. 

XX Integer counter between 0 to 5. Represents the number of keyed in 

characters while inputting set points into the controller. 

ZADIN Analog to Digital In software command. 

ZASE. ADDRESS Sets the DT2801 series board base address. 

ZATA.REGISTER Sets the data in register and data out register. 

ZCLEAR Digital to analog software clear command. 

ZDCHNL Analog to digital conversion channel number. 

ZDGAIN Analog to digital conversion gain switch at software level. 

ZEMP Data register output for digital to analog conversion. 

ZIGH High byte of analog to digital conversion output. 

ZOL# Analog to digital conversion output normalized to 5 volts. 

ZOLT# Analog to digital conversion output as a 16 bit binary number. 

ZOMMAND.REGISTER Register for the analog to digital commands. 

ZOW Low byte of analog to digital conversion output. 

ZSTOP Defines the hex value for a stop command, ZSTOP. 

ZTATUS.REGISTER Register for the analog to digital status. 

&H1 Hex value for ZCLEAR. 

&H2 Hex value for WRITE.WAIT. 
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&H5 

&HC 

&HF 

&H2EC 

&H224 

&H225 


Hex value for READ.WAIT 
Hex value for ZADIN. 

Hex value for ZSTOP. 

Base address on the PC/AT microcomputer for analog to digital conversion. 
Port number on the PC/AT microcomputer for digital to analog conversion. 
Port number on the PC/AT microcomputer for digital to analog conversion. 
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